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Electrons subjected to a strong spin-orbit coupling in two spatial dimensions could form fractional 
incompressible quantum liquids without violating the time-reversal symmetry. Here we construct a 
Lagrangian description of such fractional topological insulators by combining the available exper- 
imental information on potential host materials and the fundamental principles of quantum field 
theory. This Lagrangian is a Landau- Ginzburg theory of spinor fields, enhanced by a topological 
term that implements a state-dependent fractional statistics of excitations whenever both particles 
and vortices are incompressible. The spin-orbit coupling is captured by an external static SU(2) 
gauge field. The presence of spin conservation or emergent U(l) symmetries would reduce the 
topological term to the Chern-Simons effective theory tailored to the ensuing quantum Hall state. 
However, the Rashba spin-orbit coupling in solid-state materials does not conserve spin. We pre- 
dict that it can nevertheless produce incompressible quantum liquids with topological order but 
without a quantized Hall conductivity. We discuss two examples of such liquids whose description 
requires a generalization of the Chern-Simons theory. One is an Abelian Laughlin-like state, while 
the other has a new kind of non-Abelian many-body entanglement. Their quasiparticles exhibit 
fractional spin-dependent exchange statistics, and have fractional quantum numbers derived from 
the electron's charge and spin according to their transformations under time-reversal. In addition 
to conventional phases of matter, the proposed topological Lagrangian can capture a broad class of 
hierarchical Abelian and non-Abelian topological states, involving particles with arbitrary spin or 
general emergent SU(N) charges. 
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I. INTRODUCTION 

Quantum Hall effects are the best known experimen- 
tally observed manifestations of electron fractionalization 
above one spatial dimension^— . It is believed that sim- 
ilar fractionalization is also possible in materials with 
strong spin-orbit coupling that realize a new class of 
topological insulators (TIs) with time-reversal (TR) sym- 
metry. All currently known TIs are uncorrelated band- 
insulators^—. Quantum wells made from these mate- 
rials feature electron dynamics that somewhat resembles 
integer quantum Hall states, most notably by exhibiting 
protected gapless edge modes^. However, in addition 
to respecting the TR symmetry the new TIs differ from 
quantum Hall systems by the character of their spectra 
and by lacking the conservation of "charge" whose role is 
played by the electron's spin. The latter prevents observ- 
ing a quantized spin Hall conductivity in two spatial di- 
mensions and reduces the number of stable uncorrelated 
insulating quantum phases from infinity to only twoi^. 

The subject of this paper are strongly correlated TIs in 
two spatial dimensions whose excitations carry a fraction 
of electron's charge and exhibit unconventional exchange 
statistics^—. This research is motivated both by the 
fundamental quest for unconventional quantum states of 
matter and by potential future applications in spintron- 
ics and quantum computing 2 —. The central problem we 
address is the classification of topological orders in the 
ground states of interacting particles. We loosely define 
topological orders as distinct manifestations of macro- 
scopic many-body quantum entanglement that cannot 
be altered by tuning topologically unbiased Hamiltoni- 
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ans without going through a quantum phase transition. 
Topological orders are expressed in the phenomena such 
as the fractional exchange statistics of quasiparticles and 
the ground-state degeneracy without symmetry break- 
ing on non-simply connected spaces. We will argue that 
novel kinds of topological order are made possible by the 
Rashba spin-orbit coupling in TIs, whose description re- 
quires a generalization of the Chern-Simons (CS) effec- 
tive theory. Our main goal is then to construct a more 
general topological field theory that can capture a suf- 
ficiently broad spectrum of conventional and topological 
orders. We will discuss examples of spin entanglement 
that have no analogue in fractional quantum Hall states 
(FQHS), but a systematic classification of such states is 
beyond the scope of this paper. 

Fractional TIs can exist in various systems, and likely 
will be observed in the foreseeable future. There are 
at least three prominent approaches to obtaining frac- 
tional TIs in solid state materials. The earliest one 
relies on Coulomb interactions to facilitate spin-charge 
separation in materials with geometrically frustrated lo- 
cal magnetic moments 2 ^—. Electrons can be fraction- 
alized into neutral spinons and spinless charge-modes 
without a spin-orbit coupling, but gapped spinons can 
additionally exhibit the TI dynamics in the presence of 
a spin-orbit coupling. A more recent approach explores 
lattice models with fractional excitations, the so called 
Chern insulators^—. Such models can be TR-invariant 
and generally rely on narrow bands in the electron spec- 
trum to create favorable conditions for fractional ground 
states. A few proposals of materials that could realize 
these models have been mad o 33 i 34 . 

The third approach is to use the currently available 
band-insulating TI materials and artificially induce elec- 
tron correlations by a proximity effect in a hctcrostruc- 
ture device. For example, a conventional superconductor 
placed in contact with a TI quantum well can induce su- 
perconductivity or leave behind an insulating state inside 
the TI. A superconductor-insulator quantum phase tran- 
sition inside the TI can be tuned by a gate voltage, and 
it turns out that it would belong to the bosonic mean- 
field or XY universality class in the absence of the spin- 
orbit coupling 3 ^. This quantum critical point is sensitive 
to perturbations, and correlated "pseudogap" topologi- 
cal states can be born out of its quantum critical fan as 
a result of the spin-orbit coupling 36 . Candidate states 
are fractional TIs of spinful p-wave Cooper pairs whose 
existence is allowed by the TIs orbital degrees of freedom 
and low-energy dynamics enhanced by the spin-orbit cou- 
pling. A similar correlated TI of excitons could be envi- 
sioned in the device proposed by Seradjeh, et al£l. 

Another promising system are ultra-cold gases of 
bosonic atoms trapped in quasi 2D optical lattices. Su- 
perfluid to Mott insulator transitions can be easily ar- 
ranged to remove any energy scales that could compete 
with the spin-orbit coupling 3 ®, and thus create similar 
conditions as in the proximity effect heterostructures. At 
the same time, the recent development of artificial gauge 



fields for neutral atoms, created by stimulated Raman 
transitions between internal atomic states, has not only 
introduced the effective spin-orbit couplings 3 ^, but also 
looks very promising for generating locally enhanced flux 
densities needed for fractional state a 40 i 41 . 

The current theoretical studies of two-dimensional 
strongly correlated TIs are based entirely on adapt- 
ing the well-known descriptions of FQHS to the TR- 
symmetr y 15 i 20 i 30 ' 42 ~ — . This approach is certainly well 
motivated, but the fact is that no experimental obser- 
vations of fractional TR-invariant TIs have been made 
to date. Quantum Hall systems are sufficiently different 
from the spin-orbit-coupled materials that we must ques- 
tion their validity as an experimental basis for the com- 
plete theory of fractional TIs. Specifically, we will argue 
in this paper that the Dirac spectra of surface electrons 
in TIs pave the way to topological orders that cannot be 
fully captured by the standard CS effective theory, which 
is better suited to systems with Landau levels. 

We will instead view the two-dimensional TIs as man- 
ifestations of the SU(2) "quantum Hall physics", cre- 
ated by an SU(2) "magnetic field" that implements the 
spin-orbit coupling 4 ^. Some idealized SU(2) incompress- 
ible quantum liquids are quantum Hall states because 
they exhibit a quantized Hall conductivity of spin cur- 
rents. The unavoidable spin non-conserving perturba- 
tions in materials, however, destroy this spin-Hall effect 
and leave behind just a Z2 topological invariant when 
the TR symmetry is preserved. Interestingly, the non- 
commutative character of the SU(2) gauge fields enables 
incompressible quantum liquids without a quantum Hall 
effect even in the absence of such perturbations. All 
presently available two-dimensional TIs can be viewed 
as the non-quantum-Hall analogues of "integer" quan- 
tum spin-Hall states, where the Rashba spin-orbit SU(2) 
"magnetic flux" creates a Dirac rather than a Landau- 
like electron spectrum. The main purpose of this paper 
is then to construct and begin exploring a new topo- 
logical quantum field theory that can naturally describe 
both the fractional quantum-Hall and non-quantum-Hall 
states. It will help us predict the topological orders which 
may be specific to the Rashba spin-orbit coupling. We 
will take guidance from the experimentally established 
facts about the available TI materials and construct a 
theory that can address all of the above candidate sys- 
tems for TR-invariant fractional incompressible quantum 
liquids. 

This topological field theory will have a general form 
applicable to interacting elementary particles with arbi- 
trary charge and spin, whose dynamics is restricted to 
two spatial dimensions and affected by any type of spin- 
orbit coupling or magnetic field or both. We will use it 
to show that correlated TIs can feature excitations with 
fractional charge, spin and exchange statistics, despite 
the spin non-conservation. We will make predictions 
about the fractional excitation quantum numbers in rela- 
tion to symmetries, possible symmetry breaking, as well 
as Cooper or exciton pairing in topologically-enhanced 
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ground states. We will demonstrate the relationship of 
this topological field theory to the standard CS gauge 
theories, and point to limitations of the latter to ade- 
quately model all possible TR-invariant TIs. For the pur- 
pose of focusing on bulk topological orders, we will view 
all states of interest here as fractional TIs regardless of 
whether they have gapless edge states or not, and thus 
depart from the terminology introduced in Ref^ 5 . We 
will set up the formalism for analyzing both the Abelian 
and non- Abelian topological orders, including hierarchi- 
cal states and incompressible quantum liquids specific to 
SU(iV) fluxes and the Rashba spin-orbit coupling, which 
may have no analogue in the quantum Hall states. The 
field theory we propose can also describe conventional 
quantum phases, and possibly the universal aspects of 
phase transitions to topological states. It may be able 
to provide a broad classification scheme for topological 
states of quantum matter, analogous to that provided by 
Landau-Ginzburg theories of symmetry-broken states. 



A. Preliminaries 

Understanding complex emergent phenomena directly 
from microscopic models can be extremely difficult. It 
is often much more practical to study emergent phenom- 
ena in effective theories that specialize to the low energy 
parts of spectra. An effective Lagrangian can be con- 
structed in the continuum limit by introducing field op- 
erators to quantize the classical equations of motion, and 
by collecting all combinations of fields that respect the 
required symmetries. This method, pioneered in high- 
energy physics and the theory of critical phenomena, is 
the basis of the present analysis. This section describes 
the effective field theory of correlated two-dimensional 
TIs that we propose, outlines the principles of its con- 
struction, and surveys the structure of the paper. 

The proposed theory will be written in several different 
forms throughout the paper, but the initial discussion 
will be based on the following imaginary-time Lagrangian 
density C = £ L g + Ct'- 
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We will use Greek indices fx, v, A G {0, x, y} for space-time 
directions, Latin indices k £ {x, y} for only spatial di- 
rections, and Einstein's notation for sums over repeated 
indices. The fields ip are complex spinors with 25+1 
components whose relationship to physical spin-5 parti- 
cles is established by a duality mapping. Therefore, the 
"matter fields" ip in this Lagrangian represent vortices 
of the physical particle currents. The Landau-Ginzburg 



part £lg is the continuum limit of a standard dual theory 
of lattice bosonic particles^Zr— , adapted to the presence 
of internal (spin) degrees of freedom. Densities and cur- 
rents of particles are represented by the temporal and 
spatial components respectively of the flux $g M associ- 
ated with the gauge field matrices B^. There are 25 + 1 
independent modes of particle fluctuations that corre- 
spond to different states of spin projection on some axis 
and define the basis vectors for the matrix representa- 
tion of B^. The Lagrangian is, however, written in the 
representation-independent form. The dynamics of par- 
ticles is governed by the Maxwell term Cm in this theory. 
If the particle spin were conserved, the Maxwell term 
would have the standard non-compact form: 



Cm — 
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where Q is a coupling matrix. However, the realistic spin 
non-conservation in materials requires that certain com- 
binations of B^ modes have compact dynamics. In either 
case, particle charge and spin densities are allowed to 
fluctuate near the average values specified by the matrix 
0o, while the average current densities are zero (0; = 0). 

The topological term Ct, allowed by symmetries, 
shapes the quantum kinematics of dual topological de- 
fects in the tp field configurations. It is inconsequential in 
conventional phases such as superconductors and Mott- 
insulators, but affects the quasiparticle statistics in in- 
compressible quantum liquids. The static U(l)xSU(2) 
gauge field A^ implements any combination of external 
electromagnetic fields and spin-orbit couplings. Its com- 
ponents are SU(2) matrices, A^ = + A a ^"f a , where a M 
and A° are scalars and j a are three SU(2) generators 
in the spin-5 representation (angular momentum matri- 
ces; a € {x,y,z}). The flux components of non- Abelian 
gauge fields are 



e^ x (d v A x -iA v Ax) 



(3) 



in the matrix representation. The temporal "magnetic" 
component $o of the external flux density is inserted in 
the topological term Ct to ensure its adequate trans- 
formation under TR is the Levi-Civita tensor in 
(2+l)D space-time). There are four possible insertion 
points, and C t is symmetrized with respect to them us- 
ing anticommutators (braces). If all components of A^ 
commute with each other, then the topological term can 
be reduced to the CS form when the phase fluctuations 
in the spinor ip components drive the dynamics. It turns 
out, however, that the A^ appropriate for solid state TIs 
have non-commuting components. 

The topological term Ct is the main new ingredient in 
a field theory of this kind and we will devote most of 
the discussion in this paper to its derivation and conse- 
quences. We will derive it using the same field-theoretical 
principles that yield the CS theories of FQHS, but ap- 
plied in the context of spinor rather than gauge fields. 
The standard effective field theory of FQHS is a pure 
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gauge theory in which the CS coupling acts as a topologi- 
cal term that implements a fractional exchange statistics. 
The CS theory is constructed from the requirement that 
the action be stationary when the Hall conductivity and 
incompressible electron density are quantized as observed 
in FQHS experiments 51 . This requirement can be stated 
in an alternative form. Electrons in mutually perpen- 
dicular electric E and magnetic B fields generally have 
classical cyclotron trajectories whose orbit centers move 
at the constant velocity v = E x B/|B| 2 . The result- 
ing drift current is precisely reproduced by the kinematic 
equations of motion that make the CS action stationary. 

We will seek the analogous drift currents of spin- 
orbit-coupled particles in the section III Al Our starting 
point will be the minimal model Hamiltonian of two- 
dimensional topological band-insulators, which is by now 
well established experimentally. This Hamiltonian can 
be written in the form that couples electrons to an exter- 
nal static SU(2) gauge field with a finite "magnetic" flux 
density. We will derive the time-evolution of various cur- 
rent operators in the Heisenberg picture from a generic 
Hamiltonian of this type. The obtained equations of mo- 
tion have a direct classical interpretation according to 
the Ehrenfest's theorem. We will extract from them the 
topologically protected features of dynamics in the com- 
bined U(l)xSU(2) "electric" and "magnetic" fields. We 
will discover that topologically quantized constant drift 
currents and Hall effects are possible only when the ap- 
propriate gauge charges (charge and spin) are conserved. 
Focusing first on this special case, we will show in the sec- 
tion III Bl that the drift component of motion agrees with 
the stationary action condition applied to the topologi- 
cal term of ([T]). This will justify Ct as an effective field 
theory of quantum Hall and spin-Hall effects that can re- 
place the CS theory (later in the paper we will separately 
show how Ct can describe fractionalization, hierarchical 
quantum Hall states, etc.). 

If the above had been our only goal, we would have 
been able to construct the effective topological La- 
grangian Ct by directly considering the spin-Hall conduc- 
tivity. However, we wish to also describe the topological 
states that feature no quantum Hall effect, and only the 
SU(2) symmetry can guide us. The drift current analysis 
helps us to transparently construct Ct for any representa- 
tion of any gauge symmetry group that allows a quantum 
Hall effect. Having spinors, we can easily describe par- 
ticles of arbitrary charge and spin moving in any com- 
bination of U(l) electromagnetic, SU(2) spin-orbit and 
other fields. Now that Ct formally has the SU(2) gauge 
symmetry in any desired representation, we can directly 
apply it to the fractional TIs with the non-commutative 
gauge fields of the Rashba spin-orbit coupling. The en- 
suing spin non-conservation ruins the quantum spin-Hall 
effect, but any incompressible quantum liquid in which 
individual particles become microscopic cyclotron vor- 
tices will have excitations whose fractional statistics is 
topologically protected and correctly captured by C t . We 
will, therefore, have a tool that is more general than the 



CS theory and capable of handling the Rashba spin-orbit 
coupling. 

The construction of Ct will make it clear that the fields 
ip are not the ordinary field operators of particles, but 
rather the dual field operators that represent vortices. 
Their dynamics is provided by the non-topological part 
£lg 01 © m the most general Landau-Ginzburg form 
allowed by symmetries. The gauge field B implements 
the Magnus force on vortices in this language. Note that 
Clg governs the dynamics of smooth configurations of ip, 
while Ct is sensitive only to singular configurations. For 
this reason, the stationary action conditions for Clg and 
C t are essentially independent. The topological theory 
dual to dU) and expressed in terms of the particle field 
operators is given by the Lagrangian density ([71]) . 

The quantization of classical equations of motion, 
which we apply to obtain the quantum field theory of TIs, 
rests upon knowing the exchange statistics of elementary 
objects. The Lagrangians ([T]) with and without the topo- 
logical term Ct quantize the same classical system using 
different state-dependent exchange statistics. Therefore, 
the role of Ct is to specify quantum statistics and other 
fundamentally non-classical aspects of dynamics. This is 
done in a manner that depends on the presence and den- 
sity of topological defects in the field i\) configuration, 
which is a desirable property of a general theory that 
should describe fractional ground states. In contrast, the 
CS theory has a rigid implementation of statistics, spe- 
cific to only one particular fractional ground state. 

The technical part of this paper begins with an intro- 
duction to the SU(2) gauge-field description of spin-orbit 
couplings in the section III Al A simple single-particle 
quantum mechanics is used there to establish the equa- 
tion of motion for electrons in external electromagnetic 
and spin-orbit fields. The following section Hi Bl explains 
how symmetries and the equations of motion can be used 
to construct the field theory ([T]), and especially its topo- 
logical term. 

The next major subject of the paper are certain es- 
sential properties of the theory and its initial predictions 
in the context of the simplest Laughlin-type topological 
orders. We will first apply fundamental principles in the 
section IIII Al to show that fractionalization is mandated 
in a class of correlated quantum states that bridge be- 
tween the phases of maximally localized and maximally 
delocalized particles. These include the fractional TIs 
without spin-conservation that could arise in solid-state 
materials. Then, we will discuss in the section Ull Bl how 
and in what special circumstances the CS theories arise 
from (JTJ as effective descriptions of fractionalized states. 
Such circumstances are not met in the currently available 
TIs, and we will identify in the section Ull CI a special dy- 
namical symmetry of the Rashba spin-orbit coupling that 
can lead to new but utterly fragile topological quantum 
phases beyond the pure CS description. The following 
section IIII Dl verifies the existence of topological order 
in all these phases by calculating their ground-state de- 
generacy on a torus and other non-simply connected sur- 
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faces. Finally, the section [HIE] formally derives the topo- 
logical field theory of physical particles dual to ([T]), and 
takes a bigger perspective on the relationship between 
quantum Hall states and conventional phases of matter. 
We will touch upon the possibility of revealing the ori- 
gins of fractionalization in dynamics. The stability of 
topological orders against perturbations that violate the 
SU(2) gauge structure is briefly discussed in the section 
IIIIFI from the duality point of view. 

The following segment of the paper goes beyond the 
Laughlin-type topological orders and explores the ways 
in which the proposed effective theory (JTJ) can be general- 
ized to describe arbitrary Abelian hierarchical quantum 
Hall states (section IIV Aj) and many non- Abelian ones 
(sections IIV Bl and IIV C|) . We will demonstrate how the 
generalized topological term of (JTJ) can shape unconven- 
tional quantum statistics once the dynamics governed by 
the Landau-Ginzburg part selects appropriate low-energy 
fluctuations of the spinor fields. We will discuss in greater 
detail a class of novel and robust non- Abelian topologi- 
cal orders that can arise specifically due to the Rashba 
spin-orbit coupling. All conclusions and an outlook of 
the many remaining issues are summarized in the section 

Some aspects of this paper are variations on the well 
established results in literature, most notably the deriva- 
tions of duality and topological ground-state degener- 
acy. The theory ([T]) is closely related to the CS Landau- 
Ginzburg Lagrangian of Wen and Niu from the Ref^, 
but different enough to deserve new scrutiny. Our goal is 
both to demonstrate that (fT]) is a proper generalization 
of the CS theories, and to predict from it the possible 
new topological orders of particles with arbitrary inter- 
nal degrees of freedom. 

II. EFFECTIVE THEORY OF TOPOLOGICAL 
INSULATORS 

A. Classical and quantum mechanics 

The simplest model of a quantum well made from a 
topological band-insulator material such as Bi2Se3 or 
Bi2Te3 is given by the Hamiltonian: 

ff = «z(Sxp)r 2 +AT I -/i. (4) 

The four-component spinor wavefunction ip(r) captures 
electron's internal states labeled by the spin projection 
a z , and the orbital index t z that can be interpreted as 
the top or bottom surface of the quantum well. The vec- 
tor spin operator is S — ^a a r a , a g {x, y, z}, where a a 
and r a are Pauli matrices that act on spin and orbital 
degrees of freedom respectively (we set h = 1). In a bulk 
crystal, the two surfaces would be far apart and decou- 
pled (A = 0), so their energy spectrum E(p) = ±^vp — fi 
would contain massless Dirac states with the helical cor- 
relation S = ±ip x z between momentum and spin. As- 
suming that the chemical potential n were placed well 



within the bulk bandgap, the above Hamiltonian would 
then consistently describe the low-energy part of the full 
spectrum that contains only the surface states. However, 
electrons in a quantum well can tunnel between the two 
surfaces (A ^ 0), which opens up a gap in the Dirac spec- 
trum of surface states. A two-dimensional band-insulator 
can be obtained by pushing fi into this tunneling bandgap 
in a gated heterostructure. The given Hamiltonian is 
the minimal model of 2D electrons that both experience 
a spin-orbit coupling and have a finite bandgap with- 
out violating the TR symmetry. It has identical spec- 
trum to the model of HgTe quantum wells introduced by 
Bernevig, et al££, and may be considered different from 
it only by the choice of representation. Experimental evi- 
dence for the validity of this model comes both from bulk 
systems and quantum well o 13 ' 54 " — . 

The Hamiltonian (j4]) is related to a gauge theory for 
electrons in a static external SU(2) gauge field A. Con- 
sider: 

H> = (P-^l! + At* ,/ (5) 
2m 

where 

A = -mv(z x S) , g = t z . (6) 

The SU(2) charge g operates in the orbital subspace, 
and A^ are SU(2) matrices derived from spin opera- 
tors. Gauge transformations are specified by three angles 
9 a (r,t) combined into an SU(2) transformation matrix 
W(r,t): 

W = e lB ^ a , -> WA„W^ + -Wd^ . (7) 

For spin S = | particles, the SU(2) generators are 
7° = \o a ■ The Hamiltonians (UJ) and (|5|) produce the 
same operator equation of motion dj/dt = i [-£/", j] in the 
Heisenberg picture for the current (velocity) operator 
i=i[H,r}: 

^ = i[H, ji] = \v 2 e ijPj a z - Ave^a^y . (8) 

The symbol = eoy is the 2D antisymmetric tensor 
that implements vector cross products in the Einstein 
notation. This equation of motion (written at t = 0) il- 
lustrates the quantum cyclotron dynamics of electrons in 
its dependence on the spin and orbital index r z , which 
is the fundamental origin of all topological properties. 
In particular, we can immediately see the tendency of 
acceleration dj/di to be perpendicular to the particle's 
momentum p, and its dependence on the spin a z that 
embodies the TR-invariance. Note, however, that spin 
precession is not properly taken into account here (will be 
in the subsequent analysis). From the gauge theory per- 
spective, the cyclotron dynamics is caused by the pres- 
ence of a finite SU(2) "magnetic" flux density: 

= e^iduAx - igA v A x ) = \{mv)H^ r z a z . (9) 
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The non-Abelian nature of SU(2) gauge fields allows a 
finite flux even when the gauge field is uniform. Gauge 
transformations (0 merely rotate the flux in a spatially 
dependent way, — > W^W^ . 

The model §5§ is different from Q by the extra p 2 /2m 
term and a constant. The mass parameter to deter- 
mines the curvature of electron band-dispersions E(p) 
at larger momenta, which is indeed seen in ARPES 
experiment s 54 ! 55 . Therefore, we can regard ([5]) as a more 
accurate description of realistic systems than and 
take advantage of having the parameter to to define the 
cyclotron frequency scale and flux density. This will 
prove extremely useful in building the topological field 
theory of correlated TIs. We must note, however, that 
realistic systems do not have the SU(2) gauge symmetry. 
Still, their topological properties can be protected as long 
as the perturbations to H' that violate the gauge sym- 
metry do not remove the SU(2) flux. We shall postpone 
the discussion of gauge-symmetry violations to the sec- 
tion IIII Fl and focus first on the pure charge and spin Hall 
effects. We will explore the combined effects of spin-orbit 
couplings and external electromagnetic fields on any par- 
ticles by generalizing the Hamiltonian to the U(l) x SU(2) 
symmetry group with arbitrary spin S representation. 

Band-insulating solid state TIs exhibit a particular re- 
alization of an SU(2) "magnetic" field. We will consider 
more general situations in the following, described by the 
Hamiltonian of particles that have both an electromag- 
netic U(l) charge e and spin-orbit SU(2) charge g: 

^o= (P ' ea 2m gAY)2 ^ a o- g Agf. (10) 

This is sufficient for analyzing the cyclotron motion that 
stands behind all topological phenomena. We will im- 
plicitly assume the existence of internal degrees of free- 
dom and microscopic features that are necessary to open 
a topological gap and stabilize a TI ground state. The 
general U(l)xSU(2) gauge field A^ — a M + A a ^ a carries 
flux: 

<F = e^ x (d v A x - igAvAx) = ^ + $^ 7 a , (11) 

Its U(l) and SU(2) parts will be labeled by lowercase and 
uppercase symbols respectively, and 7° for a g {x, y, z} 
are the three SU(2) generators (angular momentum op- 
erators) in any spin-S* representation. The temporal $° 
and spatial flux matrices correspond to "magnetic" 
B and 90°-rotated "electric" Ei fields respectively, which 
together form the field tensor F^ v : 

B = F xy — —F yx , Ei = Foi = —Fio 
F[iv d^Ai/ di/A^i Ay\ — f/ii/A^* 

$m = l e ^F vX . (12) 

The traces of $ M contain the U(l) electromagnetic fields 
</> M , while their traceless parts contain the analogous spin- 
dependent SU(2) fields. Only the eigenvalues of $ M (and 



Fimv) are gauge-invariant. Defining the charge j M and spin 
J° current operators, 

jo = 1 , ji = —(Pi - eat - gAfj a ) 
to 

J$=J a , J? = \{l a ,Ji} (13) 

we obtain the following Heisenberg equation of motion 
for the spatial current components from (1101) : 

^ = i[H 0) ji] = -e ivX jAx + ^iux{3 v ,ni a } ■ (14) 
dr to Im 

One should keep in mind that all operators in this equa- 
tion are expressed at time t in the Heisenberg picture, in- 
cluding the flux operators -> e lH ot$^ e -iHo< w hich 
will process if Ho and $ M do not commute. We will 
solve this differential equation for ji(t) treated as a ma- 
trix function of time. The expectation value (ip\ji(t)\ijj) 
calculated from the solution ji(t) in any state \ip) will 
properly reflect the quantum time-evolution of currents, 
as well as the behavior of an equivalent classical system 
according to the Ehrenfest's theorem. 

Let us first consider the special case of spin-conserving 
gauge fields A^ = cl^+A z 1 ^ z whose components commute 
with each other. The resulting flux operators commute 
with the Hamiltonian, so that $ M (t) — $ M = const in the 
Heisenberg picture. Writing 

]% {t) = \°S) + K{t)i a (is) 

and organizing the scalars A" and A" into an eight- 
component vector X(t) reduces to linear differential 
equations with constant coefficients whose matrix form 
and solution are: 

^ + A\ = b (16) 

\(t) = e~ Ai A(0) + A~ 1 b . 

All eigenvalues of the matrix A are purely imaginary 
and thus generate cyclotron oscillations. The resulting 
Heisenberg current operator is 

j i (t) = c i e i » t e i ~'* u " t + 6j i , (17) 

where the first term describes cyclotron motion with 
frequencies u) ~ e<po/m, ui z — g&o/m and amplitudes 
c y = ic x appropriate for circular classical trajectories. 
The second term Sji is a constant drift current perpen- 
dicular to both "electric" and "magnetic" fields. Note 
that 5ji is state-independent and thus topologically pro- 
tected, unlike the cyclotron orbit amplitudes Cj. 

We will now concentrate on the drift current kinemat- 
ics. Setting dji/dt = in (fTH) and c ?; = in (fTT)) . we 
easily find: 

_ 1 2S 

(18) 
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It is not hard to recognize that this equation indirectly 
describes the quantum Hall effect. The amount of drift 
current is completely determined by the "magnetic field" 
(perpendicular to the sample's plane) and the in-plane 
"electric field" perpendicular to the current flow. The co- 
efficients m,k can be calculated by expanding both sides 
of this equation in the powers of 7 Z , and noting that there 
are only 25 + 1 independent matrices among (7 Z )" in the 
spin S representation. We will not pursue this expansion. 
Instead, we will need a slightly different formula 



(19) 



which is obtained by inserting e(j) +g<f>Q'y z and its inverse 
into 1(15)1. 

Now let us briefly consider the analogous dynamics 
of Rashba spin-orbit-coupled electrons. The gauge field 
© produces the Hamiltonian (ITU1) that does not com- 
mute with the flux operators ©. Consequently, the 
Heisenberg-picture operator 3>^(i) in (fLT| has a non- 
trivial time dependence. The proper way to evolve the 
gauge field operators is to treat the time evolution in the 
Heisenberg picture as a generalized gauge transformation 
that leaves all equations of motion invariant and ensures 
that all measurable (gauge-invariant) observables evolve 
according to 0(t) -> e iH t Oe -iH t. 



4> M = 0, but they are not topologically protected because 
they can have any amplitude independent of the fluxes. 
Therefore, this dynamics does not feature a quantum Hall 
effect. There are certain topologically-protected aspects 
of the dynamics, but they are buried in the oscillatory 
and momentum-dependent motion of electrons such as 
spin precession. 

The equivalent expressions (JTSJ) and (TT9")) are the most 
general operators that extract the topologically protected 
drift charge currents from any quantum Hall state of 
particles in the external U(l)xSU(2) "electromagnetic" 
fields. By symmetry, these expressions can be generalized 
to any SU(iV) group. The actual measurable currents of 
SU(JV) charges are state-dependent. 

The equation of motion for a hypothetical classical TI 
can be obtained from (|14p or (|2"T)| by equating the quan- 
tum expectation values of its left and right-hand sides in 
any wave-packet state. A wave-packet here has a spinor 
structure that should be interpreted as a representation 
of the classical spin orientation in some direction. The 
spin direction can precess, and the equation for that can 
be similarly derived from the time evolution of spin cur- 
rent operators in the Heisenberg picture. The classical 
trajectories generally involve spin precession coupled to 
orbital motion. 



B. Quantum field theory construction 



A^t) -> e lHat A^e- lH(>t + -e lHot d^e~ lHot (20) 

We can handle the flux precession by formally seeking 
the time-dependent operator solutions in the Schrodinger 
picture, where the flux operators are static: 



m = e lHot j'i{t)e- lHat 



(21) 



^ + i[H ,ji\ = -e ivX U x + £-e iuX {j' v , ni a } ■ 
at m 2m 

Like before, we can expand 8j[ as in (fT5)l to reduce the 
above equation to the form ((16)) . Its solution for Sj^t) 
must then be used to obtain the Heisenberg-picture op- 
erator ji(t) that properly captures the full dynamics. 

It is not useful for the purpose of this paper to cal- 
culate the detailed and complicated expression for ji(t). 
We will, however, benefit from revealing some qualitative 
properties of the dynamics shaped by the Rashba spin- 
orbit coupling. First, the residual commutator [H ,j-] 
in (|21j) introduces the momentum operator into the gen- 
eral solutions for ji(t), because the spin-orbit coupling 
is proportional to momentum but also contains the spin 
operators that do not commute with the flux ©. This 
means that all aspects of the current dynamics explic- 
itly depend on the electrons' momenta. Second, even the 
constant drift component of j[ is turned into an oscil- 
lating current in %, The only way to obtain a constant 
current that satisfies (|2"TT) is to insist on [H,ji\ =0. It 
can be easily seen that such solutions are possible when 



We now turn to interacting systems and construct a 
topological field theory that describes spin S particles 
and produces the equations of motion (| 19|) from its kine- 
matics. We will set e = g = 1 for simplicity and continue 
to rely on the full U(l)xSU(2) gauge symmetry in order 
to emphasize the essential TI physics. No microscopic 
information is available for a derivation of this field the- 
ory, so we must quantize (|19p the same way it is done in 
high energy physics. 

The Lagrangian we seek is required to respect the 
U(l)xSU(2) gauge symmetry in the continuum limit of 
current interest, as well as the translational, rotational 
(point-group) and TR symmetries unless the external 
gauge fields violate them explicitly. We wish to express 
this Lagrangian in terms of a spinor field whose in- 
ternal degrees of freedom naturally correspond to spin-iS 1 
particles. The usual approach would then be to associate 
tp' and ip with the particle's creation and annihilation 
operators respectively, and construct a second-quantized 
Lagrangian from the single-particle Hamiltonian such as 
((5]). However, this is the path to a microscopic formu- 
lation of the many-body Lagrangian in which the ele- 
mentary excitations are not fractionalized and have a 
pre-determined statistics. Extracting any emergent non- 
trivial statistics from the quantum vorticity of strongly 
interacting particles would be extremely difficult. 

Instead, our goal is to construct an effective theory 
that can capture fractionalization in quantum Hall states 
more directly. This theory must still be consistent with 
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symmetries and classical equations of motion. Being de- 
prived of the usual Lagrangian constructs, we need to 
consider topological terms that evaluate to zero when 
the field configuration is smooth (in a simply-connected 
space). The simplest one allowed by symmetries is given 
by C t in Q: 



Symmetries require that we choose: 



{d M -iA»){dv- %A V , $ }(<9a - iAx) 



+{(d„ - %Ap){d v - iA v ){d x - iAx), $0} 



if> 



(22) 



We will label the components ip s of the spinor ip by the 
spin projection s s {— S, . . . , S} on the z-axis, or the axis 
selected by the external spin-orbit flux $o- Variants of 
this expression turn out to be inadequate for our pur- 
poses. For example, omitting the gauge fields A^ would 
fail to produce the desired gauge-invariance and equa- 
tions of motion, while omitting the SU(2) flux matrix $0 
would yield undesired transformation under TR. Note 
that Ct changes sign under TR (ip s — > ip- s ) if $0 = 4>o 
contains only the U(l) magnetic field, while it remains 
invariant if <£> = 'Pq^/ z contains only the spin-orbit cou- 
pling. We use anti-commutators to symmetrize Ct with 
respect to the location of $0, and introduce an unknown 
coupling constant rj which cannot be determined from 
classical considerations. We will treat Ct alone as the 
Lagrangian that replaces the CS theory in its role to effec- 
tively describe topological orders. However, one should 
keep in mind that it cannot be a complete theory by it- 
self. It acts like a Berry's phase in the full Lagrangian 
(fTj), being imaginary in imaginary time. 

The action is stationary when the field configuration 
obeys: 



dCt 



IT] 



x (e^d^dx - i^ip + h.c. = 



(23) 



We used e^ x (d v - iA„){d x - iAx) = ^ vX d v d x ~ to 
derive this form. Clearly, the field configurations that 
satisfy 



(e^ x d v dx-i^)ip = 



(24) 



also satisfy (|2"3")1 . Note that the path-integral allows sin- 
gularities in ip for which the order of the above two deriva- 
tives matters. Only such singularities produce a finite 
contribution. For example, if ip(r) = e 10 in cylindrical 
coordinates, then -ie 0vX ip^ d v d x ip = e 0uX d v b\ = 2ir5(v), 
where 6 M =8^9 is the gauge field of a flux tube at the 
origin. Therefore, the condition applies to the ip's 
topological defects. If we are to interpret it as an equa- 
tion of motion for particle charge j p)Jl and spin J p cur- 
rents, we have no option but to express them as curls of 
certain ip field currents: 



1 — ^ vX Ft 1 

JPM — t u uJvX 
ja . rUvXp, Ja 



(25) 



Mi - 2 

J a - -i 



V-^o^VO ~ (d^)$ ip 



(26) 



Inserting the $0 factors is necessary for proper transfor- 
mations under TR: j p0 -> j p0 , j pi -> -j pil J£ -> -J p0 , 
J pi — > Jp 4 . This is a duality relationship. If (1231) are 
to represent particle currents, (|26[) must correspond to 
vortex currents. Even though the formulas (|26p do not 
transform properly under gauge transformations, we only 
care that the particle currents (1231) do. When the ex- 
ternal flux $0 is uniform and constant in time, we can 
rewrite the vortex charge current from (|26[) as: 

M = ^W^oVO - itf*o(W) . (27) 

and substitute it in (1231) to simplify the particle charge 
current: 

-ie^ip^od^dxip = -ie» vX ip^ d v d x ip ■ (28) 

This simplification comes from the fact that i/j^Qq ip 
is real and cannot expose any vortex singularities of 
ip to the double derivative curl e^ uX d v d x . Similarly, 
e^ vX {d v ip^)$a{d\ip) vanishes because its singular part re- 
duces to the sum of terms like e* J,,/X b sv b sX , where b s/1 = 
dfj,6 s are obtained from the phases 9 S of the individual 
spinor if) components (expressed in the representation 
that diagonalizes $o)- In summary: 



j PM = -ie^ x iP^ d v d x iP 



(29) 



Knowing the symmetry-restricted form of currents, we 
can interpret the equation of motion (|24[) . We empha- 
sized earlier that this is sensible only in quantum spin- 
Hall states, which require spin conservation and commut- 
ing gauge field and flux components, $° = ^ z l S az . If we 
multiply ([2"5|) from the left by ip^ (7 z )"$o for n = 0, 1 and 
extract the currents (|29|) from the obtained expressions, 
we find 



J~ 



(*5< 



^K^i + HK^ (30) 



where we introduced the symbols T n = if)* (-~f z ) n ip. These 
are the many-body particle currents in the stationary 
action state, which depend on the external magnetic 
and spin-orbit fluxes as well as the average vortex den- 
sities r„. We are now ready to show that these second- 
quantized equations of motion are equivalent to the first- 
quantized ones obtained in the previous section. We can 
interpret (|19l) as the renormalized current Sji = ji/jo 



9 



that describes a single particle 5 jo = 1. The ensuing 
many-body quantum average (ji) — (joSji) of the single- 
quantized formalism reproduces ()30|) after replacements 
((7 2 )™) r„, (jp) — > j pil . Analogous correspondence 
between equations of motion is found for all topologically 
protected drift currents, including spin currents and even 
currents with arbitrary powers of 7 Z placed in (|29[) . 

Therefore, Ct captures the kinematics of any topolog- 
ically protected drift currents in the combined U(l) and 
SU(2) "electromagnetic" fields. Conversely, j vlx and 
given by ([231 |2T)1) are only the drift components of the par- 
ticle charge and spin currents respectively. Recall that 
the CS theory is related to the classical drift motion in 
the same manner as £ t ■ The fluctuating currents are de- 
scribed by the gauge field B^ in (fT}. We will show later 
that Ct also determines the topological order of incom- 
pressible quantum liquids. Its ability to do so transcends 
the quantum Hall states that we used to derive it. 

The full many-body equation of motion (|23|) imple- 
ments the conservation of the non-drifting component of 
particle currents. However, its additional solutions be- 
yond pop exhibit spatial and temporal changes of parti- 
cle densities or currents. These are suppressed in incom- 
pressible quantum liquids by the Landau-Ginzburg part 
of (fT]). Such dynamics can be prominent only in conven- 
tional quantum phases, but then the entire topological 
term of (jTJ) is irrelevant as we will explain shortly. 



III. THE ESSENTIAL PROPERTIES OF 
LAUGHLIN STATES 

A. Fractionalization in incompressible quantum 
liquids 

The conventional quantum phases of bosonic particles 
that the Lagrangian (fT]) can describe are superconductors 
and Mott insulators. Superconductors can admit local- 
ized vortices at the expense of expelling particles from 
the nearest vicinity of vortex singularities (cores). This 
tends to marginalize the topological term C t of JT]), be- 
cause Ct thrives on having a finite vortex density ip^ip ^ 
and particle density (|29[) in the same regions of space, ac- 
cording to (|2"3"|) . The analogous conclusion holds in Mott 
insulators from the dual point of view. A Mott insulator 
is a superfluid of vortices whose smooth field ip config- 
urations cannot generate a finite Ct, except at regions 
in space where physical particles (topological defects of 
tjj) are localized. However, vortex currents are expelled 
from such regions (ip^ip 0). It takes strong quantum 
fluctuations to intermix particle and vortex densities and 
make Ct important. 

Incompressible quantum liquids are characterized by 
having an overlapping uniform particle density and dif- 
fused flux density. Both densities are incompressible and 
this can be formally stated by two conditions: (A) the 
density p a fluctuations are suppressed in all spinor com- 
ponents ip s = ^fp^ exp(i8 s ) of the vortex field ip; (B) 



vortices are locally coherent so that the phase gradients 
b st j. = 9^6 a follow the fluctuations of the gauge field B^ 
(whose spinor component curls represent dynamical par- 
ticle currents in individual spin channels) . In the case of 
bosonic particles, the first condition is realized in super- 
fluid and superconducting states, while the second con- 
dition holds in Mott insulators (vortex condensates) . By 
duality, (A) and (B) tend to be mutually exclusive. How- 
ever, quantum Hall states allow both conditions to be 
satisfied at the time and length scales that are probed in 
the following analysis. We will show that the theory (fT]) 
unavoidably gives rise to quasiparticle excitations with 
fractional amounts of electron's charge and spin when 
the ground state meets both conditions. Note that the 
condition (A) is more restrictive than necessary, since it 
leads to Abelian quantum Hall liquids of the Laughlin 
type. Generalizations to other incompressible quantum 
liquids are postponed until the section ITVl 

The condition (B) means that the "drift" currents ([23)1 
can be considered equivalent to the appropriate dynami- 
cal particle currents given by the fluxes of the gauge 
field Bfj, in (jTJ) . It allows us to use j PM from (|2"5]) to ex- 
press the amount of charge Q = j p odAo located within 
a small sample area cLAq during a very short interval of 
time. We can also define space-time oriented surface el- 
ements dAi — dlidt and express by AQ e u — jpidAi the 
amount of charge pushed through the sample's line seg- 
ment dli in the time interval dt. Now consider a quantum 
measurement of Q or AQdt- The outcome is random, 
but always equal to an integer multiple of the elemen- 
tary charge e. We can similarly extract the amount of 
spin S 

measurement. 

Let us denote by dC the oriented space-time contour 
that bounds cL4 M . The phases of the spinor ip s can have 
only integer winding numbers n s around the loop dC if 
ip is to be single- valued: 



Jp^dA^, which must appear quantized in any 



lis 



2tt 



(31) 



dC 



We can use (|25|) and (|26|) to express the measured charge 
and spin in terms of n s and vortex densities p s , which 
are kept constant by the condition (A): 



Q = j> dZ M j Vfl 
dC 



s=-S 

s 



27rn,fa + «$S)p, (32) 



dC s=-S 

Microscopic excitations are characterized by quantum 
numbers (Q,S Z ), where Q — 1 and S z € {-S,...,S} 
in units e = h = 1. The fixed densities p s surely can- 
not depend on measurement outcomes (Q,S Z ). Hence, 
we can view (|32p as a system of equations for n s , which 
are integers that depend on the measurement outcomes 
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(Q, S z ). We can solve these equations by requiring that 
only one of the 2S + 1 integers n s be non-zero for each 
microscopic excitation: 



n s (l, S z ) = m s 6 s ^ S z 



1 



2wp s ((t) + s$g) 



(33) 



The linearity of (|32p then generates the solutions 
(rt_s, . . . , ng) for general (Q, S z ) by adding the solutions 
for microscopic excitations. We see that the 25 + 1 num- 
bers m s must be integers. This imposes a restriction on 
the allowed values for densities: 



1 



Ps 



2irrriv 



(34) 



We have indeed obtained p s that are independent of 
the measurement outcomes, but depend on the external 
fluxes and a set of integers m s that must, therefore, char- 
acterize the ground state. There are no other physically 
acceptable solutions of (j32j) . 

The ground state charge and spin densities extracted 
from (l30l) are: 



s 

i P o= ■ 

s=-S 



S<f> Z n 



2irm« 



T z — 



s<t>v, 

s=~S 



s^ z 



2irm« 



(35) 

Various combinations of m s can lead to states with bro- 
ken particle-hole symmetry j p o ^ or magnetization 
Jp 7^ 0, especially in combined magnetic (f>o and spin- 
orbit $g fluxes. While these symmetries may be eas- 
ily explicitly broken independently of any quantum Hall 
physics, the above equations of state are purely a result 
of orbital motion and can describe spontaneous symme- 
try breaking when particles are relativistic and/or not 
Zeeman-coupled to external fields. The symmetry prop- 
erties of ground states are related to the quantum num- 
bers and statistics of quasiparticle excitations via m s . 

Since both charge and spin are delocalized in quan- 
tum Hall states, a finite area of the sample can contain 
any amounts of them on average. However, (|32[) still re- 
lates the amounts of charge and spin to quantized wind- 
ing numbers n s of dual- vortices. A single dual-vortex 
n s = <5 s ,<x in the spin channel a can be isolated in prin- 
ciple in some quantum state, for example the state of 
being localized inside a small area of the sample. If an 
experimentalist managed to suppress the fluctuations of 
n s in this localized state, he or she would measure on av- 
erage a fractionally quantized amount of charge SQ and 
spin 5S Z . More generally, a bundle of dual- vortices with 
arbitrary n s would look like a quasiparticle with charge 
and spin (in units e = H = 1): 



sq = y Hl , ss z = y 2L t 

* J Ul _ < J m _ 



s=-S 



(36) 



s=-S 



The fractionalization formulas (|32p and (|36p are indepen- 
dent of any fluctuations of the "quantum numbers" n s . 



A spin-orbit coupling such as © can favor quasiparti- 
cles that exist in superpositions of the above states with 
different SS Z . Further degradation of these "quantum 
numbers" can occur in the presence of perturbations be- 
yond the spin-orbit coupling that do not conserve spin. 
However, vortex excitations, which are dual to the above 
quasiparticles, can survive as protected fractional degrees 
of freedom because their "charges" are always conserved 
(see sections MI Fl and IIV Aft . 

We will show in the section UlI El that fractionalization 
is dynamically related to vortex "charge" . If low-energy 
vortices carry an integer number m of flux quanta hc/e, 
then charge fluctuations exhibit the fractionalized quan- 
tum e/m, observable for example in shot-noise trans- 
port measurements. The analysis in this section actu- 
ally exploits this fact. The fractional quasiparticles do 
not a priori have an unconventional exchange statistics 
in generic systems. However, quantum Hall states effec- 
tively bind a fractionalized amount of charge to a singly- 
quantized vortex, and the resulting quasiparticle is an 
anyon. The topological term in ([I]) regulates the statis- 
tics of these quasiparticles in the present formalism (and 
simultaneously gives rise to the ground-state degeneracy 
on a torus). Specifically, the fractional statistics is gener- 
ated by the ip^e^^d^d^dxtp part of the topological term. 
The following section will reveal that this is in fact the 
sum of CS self-couplings 



for each "gauge field" b sfl = 9 M # S in ground states with 
incompressible vortices. Assuming that b sfi follow the 
physical particle currents per condition (B), each quasi- 
particle in the spin-channel a acts as a source of 27r/m tT 
flux in the same spin channel. Consider two fractional 
quasiparticles with identical quantum numbers (SQ, 5S Z ) 
specified by the integers n s = 8 s . a . Their two-body wave- 
function £ CT (ri , T2) acquires the factor exp(i7 CT ) when they 
are exchanged, given by the statistical angle of the spin 
channel a: 



la 



IT 



(37) 



Note that m a = 1 corresponds to integer quantum Hall 
states of fermionic particles. 

For spin S = ^ particles, having no spin-orbit coupling 
$g = and choosing m±i/2 € Z such that v = 2m] 



1/2 



2m 



yields the Laughlin sequence of fractional quan- 
tum Hall states in the external magnetic field 4>q. The 
ground state particle density is j p o = v4>q/2'k, there is 
no magnetization J Z Q = 0, and the fundamental quasi- 
particle excitations carry fractional charge 8Q = v /2 and 
spin 8S Z = ±v/A. We see that spin must be fractional- 
ized just like charge, effectively reducing h by an inte- 
ger. A correlated TR-invariant TI ($q ^ 0) in zero mag- 
netic field cf>Q — exhibits the same combined spin and 



charge fractionalization when ?n+i/2 



ally, m s 



-1/2- 



Gener- 



is required if the ground state is to be 
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invariant under TR, which ties together the charge and 
spin fractions. Independent fractionalization of charge 
and spin generally requires TR symmetry breaking, even 
in the zero magnetic field. 



B. Chern-Simons theory of quantum Hall states 

In this section we derive from ([T]) a simplified effective 
theory in which the fluctuations of spinor ip amplitudes 
are neglected. The topological term turns into a CS cou- 
pling when all U(l)xSU(2) gauge fields and their flux 
matrices can be simultaneously diagonalized. This phys- 
ically corresponds to having a conserved spin projection 
in addition to the conserved charge. We will later show 
that pure CS gauge theories are not equally well suited 
for systems without this symmetry, and eventually argue 
that new topological orders could arise from the Rashba 
spin-orbit coupling in solid-state TIs. 

The spinor field ip = {ip-s, ■ ■ ■ ,ips) nas 25 + 1 com- 
plex components ip s — v / pjexp(i0 s ) in the theory of spin 
5 particles. We choose to express it in the representa- 
tion that diagonalizes the flux matrix $ = (f) + $q7 z , 
and assume for now that all gauge field matrices A^ and 
B^ are also diagonal in this representation. This will en- 
sure that all Ap, commute with each other in the given 
gauge. A diagonal B^ = diag09_s, M , ■ ■ ■ , Ps,n) is sen- 
sible only if the background particle "density" 9 M = 
diag(0_s,o, • ■ ■ i @S,o)<5m.o is also diagonal in this repre- 
sentation; then the gauge fields (3 sfl independently repre- 
sent the conserved 25 + 1 spin states of physical particles. 
When density fluctuations are small and mostly confined 
to microscopic scales, we may approximate p s = ip\ip s 
in all spin channels s by their averages, and capture dy- 
namics solely via the fluctuations of phases 6 S . In the 
continuum limit with translational symmetry, the aver- 
age p s are uniform, so the Landau-Ginzburg part of ([1]) 
becomes: 



C 



LG 



E 



s=-S 

+ 



K Ps 



(38) 



Here, q a are the diagonal matrix elements of Q in ([2"]l. 
which is the appropriate Maxwell term in the presence 
of spin conservation. The dots denote terms that depend 
only on p s , and can be used to determine p s in a saddle- 
point approximation. We will relate this to an XY model 
on a lattice in the section UlI El which couples 25 + 1 in- 
dependent XY fields 6 S to the non-compact gauge fields 
Ps/j,. The present continuum limit is best suited for de- 
scribing smooth fluctuations of 9 S , but one should keep 
in mind that the more accurate formulation of this model 
requires a lattice. 

Now we treat the topological term C t from (|2"!2)) to the 
same approximation. Ct is sensitive only to the topo- 
logical defects of the spinor field ip. Let us integrate 



by parts the left-most derivative of Ct in ([22)1 and write 
C t = C[ + 5C t , where SCt is the total derivative of a 
field bilinear. By Gauss' theorem, SCt picks monopoles 
d^ vX d v b sX ) j= Oof the "gauge fields" b Sfl = d„8 s . How- 
ever, the bulk monopole-charge density is zero, so SCt can 
contribute to the path-integral only at the system bound- 
aries. This can be seen from SCt oc d^j p ^ , where we inter- 
pret the (2+l)D divergence c^j PM — > of the conserved 
charge current as monopole-charge density according to 
(|29p. The only bulk contribution to the topological term 
comes from C[ , which is sensitive to the vortex singulari- 
ties <p sfl — e^ uX d v b sX ^ and yields CS effective theories 
in incompressible states. 

Let us organize the phase gauge fields b sfl into a diago- 
nal matrix = diag(fe SM ) whose flux is $g = diag(</>^). 
The bulk topological 6 S kinematics is captured by: 



~2 



in 



(S M -^)$ (^-^)W + ) 



(39) 



We emphasized only one of the four symmetrization 
terms in (1221) . but the other three denoted by dots gen- 
erate the same expression. 

In the 5 = representation, all spinors and matrices 
have a single component, and the external flux $o = 4>o 
can describe only the U(l) magnetic field. This immedi- 
ately yields the CS theory of a Laughlin quantum Hall 
stated, although not normalized in the standard fash- 
ion. We have seen in the section MI Al that the density 
p = po = \ip\ 2 must be quantized in quantum Hall states 
as p = (27tto0o) -1 , where m is a positive integer that 
specifies the winding number of the phase 9 associated 
to a single localized electron excitation. With this in 
mind, we may redefine the "gauge field" b^ = mc^ to as- 
sociate a single electron to one flux quantum of c^. The 
charge current (|29l) becomes 



— —f^ uX d 
~ <w v 



and 



C[ -> C C s = -Mr) 



~e^ x Cll d v c x 
Air 



(40) 



(41) 



We have assumed that the external gauge field satisfies 
e^ vX a fJ> d l/ a X — 0, which typically is the case. This is still 
different from the standard form 51 . First, the factor of 
— i is present because the Lagrangian is expressed in the 
imaginary time; converting it to real time removes —i. 
Second, there is an overall coupling constant r\ that could 
naively have any value. However, the CS self-coupling 
— T-W c m defines the quantum exchange statistics of par- 
ticles because j vll is the flux of c M . The value of Arjm 
must be an even integer to reproduce the bosonic statis- 
tics, and an odd integer to yield the fermionic statistics 
of elementary particles. Since m is already an integer, we 
conclude that r/ = j. This value has been already applied 
in (JTJ) . It should be emphasized again that j pfJj are local 
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drift currents of particles, which however are identified 
with the actual fluctuating particle currents & vX d v (3\ in 
the quantum Hall states where (|38|) is still capable of 
locking by — dyO to f3y (K is large enough). 



If particles have spin S 



then ipip' in (15!?)) is 

)], 

±|. We will here restrict ourselves to the 



a 2 x 2 matrix with elements yj p s p s ' exp[i 
where s, s' 

cases where all components of the SU(2) gauge fields and 
their fluxes can be simultaneously diagonalized. This 
implies [$0j$i] = an d allows choosing the gauge 
A) = —tijXi&j, Ai — — ^eijXj^a. Now the only non- 
diagonal matrix in (|39p is i/ji/j^ , so its off-diagonal com- 
ponents do not matter under the trace. Its diagonal ele- 
ments are equal in a quantum Hall state that respects the 
TR symmetry ($ = $§7 Z ) and we have discovered their 
quantization p s = (^m^g) -1 with positive integer m in 
the section MI Al The CS theory will have two decou- 
pled sectors corresponding to two spin states. In order 
to obtain the "background field" (BF) representation of 
the CS theor y 17 i 20 , we can decompose the "gauge field" 
matrix By — y( c Ii + c 1i aZ ) m t° the charge- like c C y and 
spin-like scalar components: 



C -> £rf = -i 



A fl V A ^p/J, fl pp /, 



(42) 

The gauge fields have been defined again to represent one 
unit of charge and spin by a single flux quantum: 



This theory predicts a quantized spin Hall conductivity 
a xy = Jpi/fii = (27TTO)" 1 , pertaining to the Laughlin 
spin-Hall liquids with the conserved S z spin projection. 

It should be pointed out that the CS theories obtained 
here are slightly different than the standard ones. The 
present CS "gauge fields" b S y are gradients of spinor 
phases 9 S , so their configurations admit only quantized 
flux loops by the requirement that exp(i0 s ) be single- 
valued. A flux quantum corresponds to the smallest frac- 
tional amount of charge or spin. In order to obtain a stan- 
dard CS theory, one must integrate out short length-scale 
fluctuations in the path-integral. The resulting coarse- 
grained CS theory can allow flux to diffuse and become 
consistent with unconstrained gauge fields. 



C. Fragile topological phases from the Rashba 
spin-orbit coupling 

The precise form of the external gauge field may de- 
termine certain features of the topological ground state 
that are not as robust as topological order, but depend 
on symmetries. For example, the SU(A^) Hall conductiv- 
ity is "topologically" quantized only if the appropriate 
SU(N) charge is conserved, while the bulk topological 
order can exist without charge conservation. Here we 
wish to shed some more light on these issues. 



CS theories describe naturally only the charge- 
conserving situations. An attempt to derive a CS the- 
ory appropriate for the gauge field © of solid-state TIs 
quickly runs into a difficulty. First note that (f5!J|) with 
fixed vortex densities p s expresses the topological term 
£ t for any external SU(2) gauge field Ay in the represen- 
tation that diagonalizes its "magnetic" flux $o- However, 
the Ay of © is not diagonal in that representation. The 
effective model for small density fluctuations contains not 
only the CS gauge fields c S y, but also explicit functions of 
9 3 which come from the off-diagonal elements of Ay . In 
other words, the resulting effective theory is not a pure 
gauge theory in this language. The phases 6 S generally 
fluctuate in a correlated state-dependent manner, so it is 
not permissible to simply neglect them or average them 
out. 

The CS theory depends on the external gauge field 
Ay only through its flux (fTTj) , which reduces to ^ = 
& vX d v A\ under the previously imposed restriction that 
all gauge field and flux components can be simultaneously 
diagonalized (and thus commute with each other) . How- 
ever, the same SU(2) flux can be obtained from different 
gauge fields that cannot be related by a gauge transfor- 
mation. Compare for example: 



Axy^gA^ , ^ vX d v A\ 

A 2 y = V®{Sy, x J V - 5y,yl X ) 



/i.O 



These two gauge fields have the same flux 



4> 



(d v A\ - igAuAx) = g$7%,o 



(44) 



(45) 



but only the first one has commuting components. The 
second gauge field is actually applicable to solid-state 
TIs. If we used A\y in the Hamiltonian 



(p-gA) 2 
2m 



(46) 



the spectrum would consist of macroscopically degener- 
ate Landau levels, while A 2 y would produce a fundamen- 
tally different Dirac particle spectrum. This indicates 
that there is no SU(2) gauge transformation that con- 
verts A\y to A 2 y, despite the fact that their fluxes are 
the same. Clearly, the eigenvalues of the flux matrices are 
not the only gauge-invariant quantities that characterize 
the SU(2) gauge fields. 

The qualitative distinction between Aiy and A 2 y can 
be related to global symmetries. The ideal SU(2) gauge 
Hamiltonian H\ with Aiy conserves the z-projection of 
spin C\ = 7 Z , while the same Hamiltonian H 2 with A 2 y 
conserves "helical spin" C 2 — p x l v — Py"f x where p x>y 
are momentum operator components. The commutators 
[Hi,Ci] = and [CijCy 7^ indicate that each Hamil- 
tonian (gauge field configuration) has a global symme- 
try, but the two symmetries are incompatible and cannot 
be established simultaneously. The two global symme- 
tries can be regarded as gauge-dependent: gauge trans- 
formations alter the Hamiltonians as Hi — > WHiW\ 
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so that the conserved operators must be transformed as 
d — > WCiW^ in order to keep all commutators intact. 
Note, however, that there is no gauge transformation that 
could convert C\ into C2, that is C\ ^ WC^W^ for all 
W = exp(i6» a 7 a ). 

Topologically ordered but otherwise featureless many- 
body ground states |0,) of the ideal second-quantized 
Hamiltonians T-L{Ai^) are bound to have the respective 
incompatible symmetries. The symmetry of Ai^ is "ge- 
ometrical", while the symmetry of Ai^ is "dynamical". 
|0i) and I O2) could also be the ground states of a glob- 
ally SU(2) symmetric Hamiltonian in different parame- 
ter regimes, in which case they would be separated by 
at least one symmetry-breaking phase transition. The 
qualitative difference between |0i) and |0 2 ) could be very 
deep, depending on the dynamics. For example, one 
would expect that the most stable topological orders 
in A\^ typically feature Abelian quasiparticle statistics, 
while A211 could dynamically prefer quantum liquids with 
non- Abelian statistics that we will discuss in the section 
IIV CI However, here we will consider the minimal pos- 
sible difference between |0i) and |0a). We will assume 
that both ground states can be characterized by the same 
quantized vortex densities p s in the s = ±5 spin chan- 
nels, according to (|34p. We will justify the validity of this 
assumption only in the section Till Fl 

The simplest minimally different ground states |0j) are 
both uncorrelated or Laughlin quantum liquids whose ex- 
citations can exhibit charge and spin fractionalization 
given by (|36l) . However, their specific symmetries af- 
fect the statistics of measurement outcomes in an ob- 
servable way. Suppose that one could prepare the sys- 
tem by exciting a particular fractional quasiparticle \qi) 
with desired quantum numbers above the ground state. 
The excitation \q\) of the ideal Hamiltonian T-L{Ax IJ ) will 
have a good quantum number S z and additional quan- 
tum number(s) derived from orbital motion. The eigen- 
states with opposite S z have the same energy but differ- 
ent orbital quantum numbers in spin-orbit-coupled TR- 
invariant TIs. Measuring S z of \q\) many times would 
produce a fractionally quantized average value (S z ) be- 
cause S z is conserved. Analogous quantum measure- 
ments of any other spin projection in \qi) with fixed or- 
bital quantum numbers would yield non-quantized aver- 
ages that smoothly depend on the orientation of the spin- 
projection axis. The macroscopically degenerate Lan- 
dau levels of W(Aiij.) offer many choices of orbital states 
in which this measurement statistics could be observed. 
The quasiparticle excitations |qr 2 ) of the ideal Hamilto- 
nian T-L{A2fj) behave differently. There, one must prepare 
1 92) hi a momentum p eigenstate, and measure the spin 
projection along the z x p axis in order to observe the 
fractional quantization of average values. All other pro- 
jection axes or orbital states would spoil the observation 
of spin fractionalization. 

A sharp distinction and a phase transition between 
the minimally different ground states 1 1 } and | O2) exists 
only when their respective symmetries are perfectly in- 




(a) (b) 



FIG. 1: Symmetry-protected edge spin-current states in the 
two SU(2) gauge fields given by (|44[) : (a) Ai^ that conserves 
S z , and (b) Aip, that conserves (z x p)S. The direction of 
spin-current flow at each edge is shown by the dashed arrow, 
and the local projection of spin that flows is indicated by 
short arrows. The shown flat and cylindrical geometries are 
different not by topology, but by whether the edge is curved 
in the plane of the system or not (z is always perpendicular 
to the sample plane). The two topological states (a) and (b) 
differ in each geometry by which spin projection flows and by 
whether the edge exerts a torque on the particle spin. Note 
that torque is perpendicular to spin in all cases, so it does 
not cause dissipation that would jeopardize the existence of 
gapless edge states. 



tact. The spin-Hall conductivity is quantized only when 
the spin projection (S z ) parallel to the external SU(2) 
flux ($ oc S z ) is conserved, so that the ideal |0i) is a 
quantum Hall state unlike | O2) - These properties are re- 
lated to the nature of edge modes, which are compared 
in the Figure [TJ Symmetry- respecting perturbations can 
alter the spin density and current flows along edges only 
in a manner that can be captured by an SU(2) gauge 
transformation that locally rotates spins. Since no gauge 
transformation can connect the edge modes and bulks of 
|0i) and IO2), these idealized ground states are different 
quantum phases. In fact, there are infinitely many differ- 
ent quantum phases of both kinds, characterized also by 
the number of gapless edge modes (spin Chern number 
in the uncorrelated lattice case of 1 1 } ) . 

Apart from this, the fractional quantum numbers and 
quasiparticle exchange statistics are the same in the min- 
imally different ground states |0j). We will also show in 
the next section that the ground-state degeneracy on a 
torus is the same in these two Laughlin quantum liquids. 
In fact, they have the same topological order because 
they are distinct only by properties that are not topolog- 
ically protected. The quantized spin-Hall conductivity of 
the ideal |0i), which qualitatively distinguishes it from 
I O2) , is easily destroyed by any perturbation that does 
not conserve S z , with or without the TR symmetry. The 
fragility of transitions between these phases is related to 
the fact that gaps can be opened in the edge state spectra 
by perturbations that remove the symmetries of the ideal 
|0i) and IO2). The TR-symmetry alone protects only a 
Z 2 edge-state distinction between topological state o 15 ' 20 . 
The lack of symmetry also spoils the bulk measurements 
of fractional spin discussed above, but does not jeopar- 
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dize the topological order expressed via the topologically 
protected numbers m s in (|34p . 



D. Topological ground state degeneracy 

Here we calculate the ground state degeneracy of frac- 
tional TIs on non-simply connected surfaces. We will fo- 
cus on the Laughlin sequence of TR-invariant fractional 
states of spin S = i particles. Our main goal is to es- 
tablish the existence of spin-related topological orders 
despite the fact that the Rashba spin-orbit coupling of 
solid-state materials spoils spin conservation. This is mo- 
tivated by our interest in topological orders that are cre- 
ated by the Rashba spin-orbit coupling, rather than the 
ones merely perturbed by it. 

Let us first review the procedure for extracting the 
topological degeneracy in the well-understood situations 
when spin and charge are conserved^. Consider the TR- 
invariant SU(2) Hall effect shaped by the external gauge 
field Ai^ in (|4"4")l . The two spin projections are decoupled 
and experience opposite magnetic fields. This admits a 
pure CS gauge theory description of low energy dynam- 
ics in quantum Hall states. The TR-invariant Laughlin 
state discussed in the section IIII Al is characterized by 
vortex densities p s = iplip s — l/(7rm$o)j where s = ±5 
for spin S = h particles, and m is a positive integer. We 
assume that vortex fields ip s — v / p7exp(i0 s ) are locally 
coherent in quantum Hall states, so their phase gradi- 
ents b sfl = d^Og are locked to the gauge field B^ that 

We may, there- 
to arrive at the 



represents the physical particle currents 
fore, integrate out S M in (JTJ and use ([39 
following effective CS theory: 



C 



cs 



E 



Airm 



Airm 



(47) 



We have defined a s — 2s = ±1, and omitted the 
generated vortex Coulomb and current-current interac- 
tions. The latter is justified by our exclusive interest 
in the ground states and the fact that vortices are al- 
ways gapped. Recall that the "gauge fields" b Sfl = d^9 s 
admit only quantized flux loops by the requirement that 
exp(i9 s ) be single- valued. This constraint enables a sharp 
distinction between gauge field configurations that corre- 
spond to the ground-state manifold and excitations. Any 
quantized flux tube is bound to cost a finite amount of 
energy due to its narrow core. The exception are flux 
tubes that reside outside of the space in which the La- 
grangian density is defined. If the space is shaped as 
a torus, a flux tube can be threaded through one of its 
openings without paying any core energy, as illustrated 
in the Figure [2ja) . 

The situation is temporarily complicated by the pres- 
ence of 6 s o 7^ in (|4"T|) . which is related to the back- 
ground density of particles. In normal circumstances, 



this nucleates flux lines that stretch along the path- 
integral's time direction at an average spatial density de- 
termined by s o. The presence of a net bulk magnetic 
flux is reflected by the corresponding circulating currents 
along the system's boundary. However, we will consider a 
torus geometry of space, which has no boundaries. The 
torus geometry frustrates the system by allowing only 
the configurations of temporal flux lines whose net flux 
is zero. A finite O s o is, therefore, fairly innocuous on a 
torus. It is better to not have any temporal flux lines 
at all then to compensate every flux line by an anti-line. 
The exception are, again, flux tubes threaded through 
the torus openings. We conclude that ground-states on 
a torus are shaped by phase 6 S configurations that wind 
an integer number of times around any torus opening. 

If the torus sizes in both x and y spatial directions are 
I, then the relevant ground-state configurations are: 



(r) = — (n sx x + n sy y) =>■ b si = 



2nn s 
I 



(48) 



where the integer winding numbers. Since b s 



dfj,9 s , we also find 
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sy I 



(49) 



and e^ vX d v b s \ = indicates that, indeed, no flux pen- 
etrates the torus surface. The configurations P5j) are, 
however, not the only relevant ones because they pro- 
hibit electric fields (spatial flux) on the torus. If the 
magnetic flux threaded through a torus opening changes 
in time, dQn si ^ 0, then electric field loops should appear 
on the torus surface according to the Faraday's law. This 
is physically required because a threaded magnetic flux 
is merely a circulating supercurrent flowing around the 
torus, so changing it requires an electric field pulse. 

The Faraday's law is a consequence of flux conservation 
in the path-integral. A flux-tube stretching in the time 
direction is a constant magnetic field, while its bending 
toward a spatial direction turns it into an electric field 
flux line that unavoidably coincides with the change of 
magnetic field. Electric field is actually perpendicular to 
the spatial flux vector, so the flux-tube bending into var- 
ious spatial directions generates circulating electric fields 
around the place where the magnetic field changes. It 
is necessary to integrate out short length- and time-scale 
fluctuations in the path-integral in order to recover the 
continuum Faraday's law from the diffusion of quantized 
flux loops. This is visualized in the Figure HJb). 

Now we can look for the missing important field config- 
urations in (|48|) . If one viewed the (2+l)D torus space- 
time as the boundary of a four-dimensional hyperspace, 
then a quantized magnetic flux tube threaded through a 
torus opening could bend in the "fourth dimension" to 
accommodate don S i ^ 0. This would cost energy when- 
ever the tube touched or punctured the torus, but the 
cost can be arbitrarily small if the entry and exit points 
are sufficiently close together. We may regard such tube 
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FIG. 2: (a) Magnetic flux loops with n x and n y quanta 
threaded through the torus openings. This is a depiction of 
the field configurations (|48p that have a vanishing energy cost 
in the limit I — > oo of infinite torus size, (b) The emergence 
of the Faraday's law in the (2+l)D path-integral. A flux line 
stretching along the imaginary time direction r represents 
the scalar magnetic field B = e 0lJ dibj (thick blue arrow). 
The magnetic field is constant in this drawing until r = to 
when it vanishes. Flux must be conserved, so the flux line 
(thin black arrow) must bend into a spatial plane in order 
for B to vanish. A spatially-stretching flux line represents 
the vector electric field Ei = di bo — dobi (thick red arrow), 
which is related to the spatial flux component by the right- 
hand-rule, e^d^b 

v — e ij Ej . The flux line may extend in any 
direction in the xy plane with equal probability. Whenever 
the path-integral dynamics leads to flux diffusion, these direc- 
tion choices are averaged so that the coarse-grained electric 
field lines form closed loops and have magnitudes fixed by the 
flux conservation to the values predicted by the Faraday's law 



-OqB. This drawing illustrates a segment of the 



torus space-time, but the analogous Faraday's law is gener- 
ated by instanton events in which the quantized flux tubes 
threaded through the torus openings in the figure (a) bend 
to change the magnetic field in time and penetrate the torus 
space. 



intrusions into the torus space-time as instanton config- 
urations of 8 S . They must be included in the low energy 
dynamics, but there is no need to construct them ex- 
plicitly. We may simply coarse-grain the action on the 
torus by integrating out the short-scale fluctuations, and 
require that the Faraday's law be dynamically obtained 
from instanton events. The coarse-graining will allow flux 
to diffuse on the torus surface, but the threaded magnetic 
flux through the torus opening remains quantized. Now, 
note that keeping b S i from (|4"5)l while setting & s o to zero 
would yield the configurations 



b s0 = , b s 



(50) 
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never renormalized). If we substitute (IBTfl) in the coarse- 
grained Lagrangian density, convert it to real time and 
integrate over the spatial coordinates of the torus, we ob- 
tain the Lagrangian of the spatially-uniform low-energy 
modes: 
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(51) 



Note that & s o drops out completely due to the torus 
geometry, and hence affects only the dynamics of exci- 
tations. The canonical coordinates are q S i = n S i, so the 
corresponding canonical momenta are: 
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(52) 



and the Hamiltonian is: 
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This is equivalent to the quantum mechanics of a single 
particle with an internal spin-state, in a spin-dependent 
external magnetic field. The particle is allowed to move 
only on a discrete square lattice whose lattice constant 
is 1. The equivalent gauge fields a S i and their magnetic 
fields &c are: 



b s = (,,(),'!, 



27T<7 S 



(54) 



The amount of flux per plaquette is ±27r/rn, so that there 
are 1/m flux quanta per plaquette. The above Hamilto- 
nian, thus, corresponds to a Hofstadter problem, whose 
spectrum is m-fold degenerate in each spin sector—. The 
total ground-state degeneracy is m 2 . 

We are now ready to analyze the main problem of inter- 
est. Consider the gauge field A-^l from (|4"4")l . which cap- 
tures the Rashba spin-orbit coupling in solid-state TIs. 
An attempt to derive the CS theory from (|39[) ends with 
the following effective theory: 
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(55) 
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that precisely implement the correct Faraday's law. 
Therefore, these are the full coarse-grained gauge field 
configurations that adequately represent the dynamics 
of the ground states. 

The coarse-grained Lagrangian density has the same 
form as (|47p . but with a renormalized coupling of the 
Maxwell term (the CS couplings are topological and 



The spin non-conservation should introduce compact 
components in the Maxwell term @, but we have ex- 
panded them to the quadratic order given that the rele- 
vant flux densities (|50|) are extremely small in the ther- 
modynamic limit I — > oo. Even though vortex current- 
current interactions are omitted again, this is not a pure 
gauge theory because S z is not conserved. Nevertheless, 
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the relevant 9 S configurations in the ground-state man- 
ifold are still given by (|48l) and instantons for the same 
reasons as before. Instantons always introduce flux lines 
into the torus space-time and hence rotate 9 S in a man- 
ner that averages exp(i9 s ) to zero over space-time. We 
have no means to capture this exactly in the complicated 
Lagrangian density (|55p . However, we can consider a hy- 
pothetical worst-case scenario in which instantons fail to 
annihilate the exp(i# s ) factors even after coarse-graining. 
If instantons and coarse-graining only produced (j50[) as 
the relevant gauge field configurations 6 S „ (which are now 
liberated from 6 S due to flux diffusion), then integrating 
out (|55[) over the spatial coordinates on the torus would 
yield the real-time Lagrangian: 



2q 2 



(56) 



(d n tx - d n ix ) + i(d n ty + d Q n iy ) 

The factors e l ^ r ~ Si " > are averaged to zero unless rifi = 
because of the windings (|38]l. Consequently, the sec- 
ond part of the Lagrangian always vanishes, so this La- 
grangian is actually identical to (|5T|) and we are tempted 
to conclude that the ground-state degeneracy on a torus 
is again m 2 . But, before we make any conclusions 
we must address the gauge-dependence of (|55|) because 
choosing a gauge different than (|48|) would invalidate the 
above averaging of the e 1 ^ -6 ^ factors. Since we are not 
dealing with a true gauge theory, different gauges are re- 
ally different physical states. Let us recall that spatially 
varying 9 S configurations generally represent vortex cur- 
rents (p?6"|) in the original theory (JXJ) . We focused ear- 
lier only on the circulating vortex currents that produce 
the quantized flux tubes of b sfi . However, any vortex 
density or current produced by a deviation from (|48j) im- 
plies the presence of vortex excitations. Vortices are fully 
gapped by the assumption that the ground state is not 
a Mott insulator (and they are gapped even in Mott in- 
sulators by the Anderson-Higgs mechanism according to 
our simplified model that neglects the physical photon 
fluctuations). Therefore, the "gauge" we used in deriv- 
ing (|56p is exactly pertinent to the ground-state manifold 
and sharply separated from other "gauge" choices by an 
energy gap. The only important deviations from (|48[) 
are instantons, but they are bound to help rather than 
hinder the averaging of e^t" ' ±> to zero. 

We conclude that the ground-state degeneracy of a 
TR-invariant SU(2) Laughlin state on a torus is always 
to 2 . The specific spin- non-conserving form of the exter- 
nal gauge field only determines the global symmetry 
of the ground state, and the details of spatial fluctuations 
in the Lagrangian density that are relevant for the exci- 
tation spectrum. We are also assuming that Aip, shapes 
the energy landscape in a way that gives rise to the spe- 
cific edge states discussed in the previous section. Quasi- 
particles are fractionalized, but their conserved spin-like 



quantum number is the spin-projection perpendicular to 
momentum according to the right-hand rule. 

A more general derivation of the topological ground- 
state degeneracy relies on the "topological symmetry" 
transformations of the effective Lagrangian in quantum 
Hall states. This procedure reveals the degeneracy on 
Riemann surfaces of arbitrary genus g. We will not dis- 
cuss it here because it essentially follows the calculation 
of the Refi^ 2 - and arrives at the analogous conclusion 
that the degeneracy of all TR-invariant Laughlin states 
is to 2s . One would have to start from the Lagrangian 
([71)1. which is dual to ([T]) and describes physical particles 
directly. Then, one would derive an effective theory of 
quantum Hall states from this Lagrangian (by fixing the 
particle densities in all spin channels in the topological 
term). This effective theory is essentially the same as 
the Landau- Ginzburg theory with a Chern- Simons cou- 
pling from the Ref.— . The only difference is in the choice 
of a (singular) gauge and an implicit constraint for the 
flux quantization of the CS gauge field, but this does 
not affect any symmetries. Therefore, the same symme- 
try analysis can then be performed to demonstrate topo- 
logical orders on Riemann surfaces. The ground-state 
degeneracy is shaped by the low-energy dynamics of par- 
ticle or vortex fields, and not directly by the specific spin 
conserving or non-conserving form of the external SU(2) 
gauge flux. 



E. Duality 

This section discusses the duality relationship between 
the vortex theory (p} and the theory that directly de- 
scribes the dynamics of physical particles (jT4")l . We 
will assume that the particles are bosons and formally 
derive duality using the XY and CS descriptions of 
dynamics^Zr— , which are not always available as we have 
seen in the previous sections. This and the symmetries 
will guide us in making a conjecture about an approxi- 
mate duality that holds generally at the level of the spinor 
Lagrangians |1| and (174")) for any type of particles. 

When the phase fluctuations of the spinor tp compo- 
nents dominate the dynamics, the Landau-Ginzburg part 
of ([TJ) can be viewed as the continuum limit of the gauged 
quantum XY model on the square lattice. Its action can 
be formulated in the discretized imaginary time: 

SvXY = --$>os(A^ s -/3 SA1 ) (57) 
Qs 

s L " r,fi 



The angles 9 S are the phases of the complex spinor com- 
ponents ip s in the representation that diagonalizes the 
gauge field matrix = diag(/3_ SiAl , . . .,Ps,n)- In this 
microscopic model, however, 9 S live on the (2+l)D space- 
time lattice sites r, and f3 sfl live on the lattice bonds, 
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where ji E {0, x, y} labels the three bond orientations 



and p s 
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crete lattice derivative is A M S)I . 
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9, r and the 



discrete curl e )J " vX Av(3 s \ is given by the oriented sum of 
PsJtt') on the bonds (rr') of a square plaquette. 

The sum of all gauge fields f3 s ^ is non-compact as im- 
plied by ((T|), because the physical particles must be able 
to exist in fully delocalized states such as supcrfluids 50 . 
Certain linear combinations of f3 s ^ are allowed to be com- 
pact since spin need not be conserved, but we will assume 
for simplicity that every /3 SM is individually non-compact 
in (|57|) . We will also neglect the electromagnetic field 
fluctuations to which electrons are coupled; otherwise, 
the sum of /3 sfl associated with charge fluctuations would 
be gapped out by the Anderson-Higgs mechanism. The 
formal duality transformation maps (l57l) into the XY 
model of particles 4 ^ whose continuum limit is given by 
the Landau- Ginzburg part of (f74|) . 

The gauged XY model ([57)l represents the dynamics of 
vortices, while its dual XY model represents the dynam- 
ics of particles. According to duality, the transition from 
the superfluid state of particles to a Mott insulator can 
be viewed as the condensation of quantized vortice o 48 ' 59 . 
Particles are mobile and coherent in the superfluid state, 
while vortices are gapped and localized into a vortex lat- 
tice if their density is finite. A Mott insulator is a dual 
reflection of the superfluid where particles and vortices 
exchange their behavior. 

We can qualitatively view incompressible quantum liq- 
uids of bosons, including quantum Hall states, as "ar- 
rested" Mott transitions in which both particles and vor- 
tices are abundant and mobile, yet uncondensed and con- 
trolled by the cyclotron scales. Duality allows simultane- 
ous mobility of both particles j pfJ , and vortices j vfl only if 
vortices are "attached" to particles and both have incom- 
pressible densities. This prevents a vortex from passing 
through a particle, and thereby prohibits strong phase 
fluctuations that would localize the particles (and vice 
versa by duality) . We must imagine that the superfluid- 
like correlations are not locally lost in an incompress- 
ible quantum liquid, but particles have nevertheless be- 
gun localizing with respect to other nearby particles. In 
the presence of external fluxes, the wavef unction of a lo- 
calized particle must acquire vorticity, and this consti- 
tutes the microscopic mechanism for "flux attachment" . 
The quantum Hall regime is achieved when every parti- 
cle becomes a microscopic "cyclotron" vortex. The mo- 
tion of particles is then governed by Lorentz and Mag- 
nus forces, whose effect is captured by the topological 
(Chern-Simons) terms in the Lagrangian. 

An example of conditions that could produce an in- 
compressible quantum liquid is found beyond a first-order 
transition from a vortex lattice superconductor to a cor- 
related insulator. Vortex lattice melting is generally first- 
order—, but this implies that the superconducting system 
enters an insulating state while its order parameter mag- 
nitude is still finite. In other words, the superconductor's 
phase coherence survives at short length-scales given by 



the separation between vortices in their liquid state. This 
is enough to define isolated vortices in the first place and 
give rise to Magnus forces on them, provided that vortices 
are small and compact as is the case when the number of 
particles per vortex is small. 

We will now include the topological term £ t of ([1} in 
the duality analysis, in order to discover its form in the 
direct theory of physical particles. Unfortunately, the 
form of Ct is precisely defined only in the continuum 
limit. This can be appreciated by considering an equiva- 
lent CS coupling e fJ,uX b ti d v b\, which couples the particle 
current j pfl cx & vX d v b\ to the vortex current j VM ~ b^. 
If a CS term were to be defined microscopically, then the 
particle current would have to live on the bonds of the 
direct (physical) lattice, while the vortex currents would 
live on the links of the dual lattice. These two lattices 
do not coincide, so the lattice CS coupling j PM j V /i is am- 
biguous. Instead of trying to deal with this issue, we will 
carry out the duality mapping in the continuum limit. 
Hence, we will give up keeping track of the detailed lat- 
tice effects on the particle and vortex dynamics, but they 
are anyway known from the exact duality mapping in the 
absence of topological terms. Our interest is mainly to 
discover the continuum limit of the topological term in 
the particle theory, and we will obtain the adaptation of 
the result from the Refi^ 2 - to the specific features of the 
present formalism. 

Our starting point is the effective theory given by ([55)1 
and (|39p. We require the presence of the U(1) 2S+1 sym- 
metry that decouples different spin channels in the di- 
agonal representation of $o- In that case, all XY fields 
9 S fluctuate independently, and we need to only focus 
on the dynamics of one of them at a time. A single 
copy of the CS term looks like Pdj) with the topolog- 
ically fixed coupling rj = where the normalization 
d^9 s = b Sjli = m s c sf i is used for the CS gauge field. 
Therefore, the effective theory of vortices in a single spin 
channel s G {— S, . . . , S} takes the form: 
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-e^ x A s ^b sX . (58) 



27rm, 



Even though b st± reflects all local fluctuations of 9 S . only 
its singular (gauge-invariant) part contributes to the CS 
coupling. It is therefore convenient to carry out a singu- 
lar gauge transformatio n 61 ' 62 and separate two kinds of 
fluctuations 9 S = 9' s + where 9' s are smooth and 9" 
are singular phase configurations: 



e^%(d x 9' s ) = , 



(59) 



We have redefined the "gauge field" b stl by fixing its 
gauge. In order for vortex spinors to be single-valued, 
exp(i# s ) must be single- valued and hence the space-time 
configurations of b sll can allow flux to be concentrated 
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only in quantized loops as if the dynamics of b sfl were ex- 
tremely compact. Monopole configurations of b sfl would 
have been allowed if only they didn't have to be com- 
pensated by p sfjl , which is non-compact. We can now 
promote b s ^ to arbitrary gauge-field configurations and 
impose the constraint of "compactness" via the Pois- 
son summation formula, by introducing an integer- valued 
fluctuating current field J sfl in the path- integral z: 

r. OO 

z = I d9' s db sfl dp sfl e s J] £ S {e^ x d„b sX - 2wn sfl ) 

sfi n s ^— — oo 

«/W/> Y>— ,60, 

An integer J sll is associated with every lattice bond in 
the microscopic lattice formulation of the path-integral, 
but here we must treat J SiL as an integer-valued vec- 
tor function of continuous space-time coordinates. The 
transformed Lagrangian (|55j) reads: 
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M + b Sfl -/3 SM ) 2 - iJ s ^" x d u b sX + ■■■ (61) 



Let us decouple the two XY-model-like terms in (|58j) 
by applying the Hubbard-Stratonovich transformation in 
the path-integral: 



exp 



dy exp ( -|- + xy ) (62) 



We need to introduce two real-valued Hubbard- 
Stratonovich current fields that appear in the trans- 
formed Lagrangian, j sll for particles and for vortices: 
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(63) 



Then, integrating out Q' s and /3 SAt multiplies the measure 
of the path-integral by the Dirac delta functions that 
implement the following constraints: 



Strictly speaking, only the transverse (flux-changing) 
modes of /3 SM are physical and fluctuating, but integrat- 



ing them out imposes the constraint j 



e^ x dj 



sX 



only on the transverse part jn, of the full vortex current 
j sfi = jjfj_+jsfi- The transverse and longitudinal currents 

are defined by restrictions d^j^ — and ^ uX d u jl x = 
respectively, so that integrating out 9' s quenches the lon- 
gitudinal current by imposing c^ji'^ — > 0. Note that 8' s 
does not include the singular configurations that produce 
a non-zero flux of 6 SM = 8^9 s , but this does not affect its 



coupling to the longitudinal current. Once the longitudi- 
nal current is suppressed, we can identify the total vortex 
current with its transverse part and solve the resulting 
constraints (|64[) by introducing a scalar and a gauge field: 
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(65) 



The scalar field <p' a is not allowed to have singular configu- 
rations because the curl of j sll must be generated entirely 
by a s>jL . Substituting the above expressions in (|63|) gives: 
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where we have introduced symbols 



(66) 
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(2irm s ) 2 Kp s 



(67) 

to simplify notation (note that n s > restricts the al- 
lowed values of Q s q and m s ). Now we are ready to inte- 
grate out bg^, the last field that originated in the vortex 
Lagrangian ([58| . This can be done by completing the 
square: 



e^ x b m d v b s x - 2e^ x b su d v 



e^ x (b s ^ - a SM ) d v (b 



l_lUl/Q>s\ 

a — a s \) 



(68) 



e^ vX a S ud v a s x 



where 
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(69) 



Integrating out 6 s/i 



a S n merely multiplies the 



path-integral by a constant and leaves behind the effec- 
tive Lagrangian of physical particles: 



+in s (d (p s 
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e^ X a slL d w a s x . (70) 



(64) Expanding the CS coupling for a yields: 



e^ x a s M sX 



e^ x a s ^a s x (71) 



- iJ s ^ vX d v a s x - nTm s e^ x J sll d v J sX . 

Since J S[1 G Z, the last term in this expression is an 
integer multiple of 2ixi for bosonic particles that require 
m s to be an even integer. Then, this term does not affect 
the path-integral and can be ignored (because e 2win = 1). 
In contrast, fermionic particles have an odd m si so for 
them the last term introduces interference of the Z 2 kind 
in the path-integral. We will not attempt to study this 
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interference here; that can be done only in the lattice 
limit, and may prove to be a formidable challenge since 
duality mappings for fermions are not known. Instead, 
we will simply focus on bosons and drop the last term in 
(|7ip . The summation over J sil has no quadratic weight 
any more and thus forces the gauge field a sfl to collect 
flux only into quantized vortex loops according to (1601) . 
From this point on, we can view a sfl as the gradient of a 
singular phase. Let us define an XY field <p s via d^ip,, = 
d^v's + a sti which is allowed to contain singularities, and 
fully eliminate a sfl from (|70p in favor of a sfl : 

+in s (d fs - A s0 ) - e^ vX a s ^d u a sX . (72) 

<j,p _ e ^ vX Q^A sX is the externally imposed static flux in 
the chosen spin channel. If we pick the gauge in which 
a sfl = d^ifs, this Lagrangian becomes entirely analogous 
to the starting combination (1551) of a continuum "XY" 
model and a CS term. Note that the first term (l/2e*) 
is redundant because it is obtained merely by taking the 
curl of the second term. Particles are here coupled only 
to a non- fluctuating external gauge field A sfl . 

The continuum limit (|72|) by itself does not define a 
charge quantum, since the field tp s can be arbitrarily 
renormalizcd. A reference for the normalization of (p s is 
provided only via the topological CS term, which spec- 
ifies the statistics of (vortex) excitations. However, the 
microscopic lattice rendition of (|72[) is compact in such 
a manner that the charge quantum is the fundamental 
electron's charge. The natural (compact) XY fields are 
the angles tp s that we wrote, judging by how they cou- 
ple to the external gauge field A sfl . If the particle density 
were an integer number per lattice site in each spin chan- 
nel, we would have n s — and the lattice XY action of 
particles without the topological term would be: 



(73) 



We are now ready to construct the continuum-limit 
Lagrangian of the physical particle spinor field £, whose 
limit of small density fluctuations is given by (|T2"|) . We 
will immediately write it for particles of arbitrary spin 
S. Defining the particle spinor £ = ■ • ■ ,£s), where 

£ s = y / 7i7exp(i(p s ), we obtain: 



(74) 



This Lagrangian exhibits the non-relativistic dynamics 
that we found in the duality mapping, but relativistic 



dynamics is also possible in special circumstances (inte- 
ger particle filling of lattice sites). The chemical potential 
matrix t determines the average particle densities n s in 
all spin channels, while the kinetic energy of particles 
in different internal states is determined by the matrix 
K = diag(i^_5, . . . , Kg) that equals a constant multiple 
of the unit-matrix in normal circumstances. The em- 
bedded factor of Oq" 1 in the topological term assures its 
proper transformation under TR. The dynamical particle 
charge and spin currents take the usual non-relativistic 
form in terms of the particle spinors: 
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(75) 



j? = -- 



It should be pointed out that the topological term 
C p t is gauge-invariant, but written in the gauge fixed by 
the non-dynamical spin-orbit coupling embedded in A^. 
Strictly speaking, the derivatives <9 M — > d^—iV^ in C pt are 
covariant, where V p is the static background gauge field 
dual to A^ and minimally coupled to the vortex spinor ip 
in (JT|). We never wrote V M before because it equals zero 
in the chosen natural gauge, but it is formally added to 
the purely transverse fluctuations of the dynamical gauge 
field Bfj, in order to carry the gauge transformations of 
Bfj, when the spinor ip is transformed. While present for 
general SU(iV) symmetry groups, the gauge-invariancc 
of Cpt can be explicitly revealed in the U(l) duality by 
merely shifting /3 sfi — > (3 sfl +v sfJj by a constant background 
gauge field u SjU , where = diag(u_s )At , . . . , 

We conjecture that the topological field theory (|74)) 
of physical particles is generally dual to the correspond- 
ing theory of vortices (JTJ), despite the fact that we could 
derive it only in the case of bosonic particles whose dy- 
namics respects the global spin \J(1) 2S+1 symmetry. The 
duality relationship between (JXJ) and (1741 is established in 
the following qualitative sense. The spin S particles are 
directly represented by the spinor field £ in (|74p . which 
is minimally coupled to the external U(l)xSU(2) gauge 
field A^ that embodies magnetic fields and spin-orbit in- 
teractions. The topological term regulates the quantum 
statistics of vortices, which are line-like topological de- 
fects in the (2+l)D space-time configurations of £. The 
dual theory (Q]) directly describes the dynamics of vor- 
tices, represented there by the spinor field ip. Each spin 
channel of particles corresponds to a spin channel of vor- 
tices. By duality, vortices must be minimally coupled 
to a gauge field B^ that represents the physical particles 
via its fluxes. The topological term of ([Tj regulates the 
statistics of physical particles. 

There are several interesting features worth pointing 
out. First of all, the lattice version of the theory ([74]) 
applied to bosonic particles can describe conventional su- 
perconducting and Mott-insulating phases. The written 
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continuum limit represents by n s = the excess den- 
sity of particles in the spin state s relative to an integer 
number per lattice siteMiSi. The dual "flux" 9^ = 9 <^, 
depends on these densities according to (|57|) and thus 
acts as a source of Magnus forces on vortices in the dual 
theory ((T|). However, the conventions used in the present 
field-theory constructions impose the same kind of sym- 
metry transformations on 9 M as on the SU(2) "magnetic" 
flux <&o. For example, a TR- invariant system will feature 
9 M = nj z S fl fl with n being the density of particles, which 
makes a distinction between the particles of opposite spin 
even when their densities are the same. The ensuing vi- 
olation of the full SU(2) symmetry is merely a reflection 
of the fact that densities are measured with respect to an 
integer filling of lattice sites. A stable ground state has 
full SU(2) symmetry only if it is a Mott-insulator with 
the same integer number of particles per site in all spin 
states, but then n = and 9o = 0. 

Magnus forces on vortices are proportional to particle 
densities rather than the dual "flux" 9o as naively por- 
trayed by (fTJ) and (|58p . This can be revealed by renormal- 
izing the gauge fields /3 SjU — > m s f3 sfM that capture ordinary 
matter fluctuations in the dual theory. Using (|67[) . the 
first two terms of ((58)) become: 

^{d^s - m s ft,) 2 + ^(£ FA 9^ S A - 27rn s <V,o) 2 • 
2 on qi 

We see that the numbers m s are vortex charges. In con- 
ventional phases m s = ±1 and there is no fractionaliza- 
tion. However, the interplay between various coupling 
constants in the Landau-Ginzburg theory (1741) can in 
principle set such densities and dispersions of particle 
modes that the charges m s become non-trivial via (|67[) . 
This is not a complete picture of how fractionalization 
arises dynamically, but it is more general than a rigid 
kinematic CS description of fractionalized states. We will 
not make any further attempts to study the dynamical 
origin of fractionalized states in this paper. 

We found in the section IIII Al that vortex densities 
p s are locked into the values of m s in quantum Hall 
states, and similar locking between the particle densities 
n s and m s can be extracted using the duality analysis 
from this section. Namely, mj 1 are filling factors, so that 
n s = <i>o s /27rm s and 9o = $o according to (1671) in quan- 
tum Hall liquids. Therefore, 9o is rigidly determined at 
least in all Laughlin quantum Hall liquids, irrespective of 
which fractional charges are selected by dynamics. The 
appearance of 9q 1 —> Q^ 1 in the topological term of (f?4")) 
is not alarming because this relationship between 9o and 
$0 does not hold in conventional states of matter that 
always become stable in the $0 — > limit. In such con- 
ventional states, 9o is determined by the particle density 
as discussed above. 

The theory (|74[) might apply to a much more general 
context than the one we derived it in. For example, we 
could formulate it for fermionic particles even though du- 
ality mapping is then not known. Furthermore, we will 
see in the sections IIV Al and IIVBI what generalizations 



are necessary for hierarchical Abelian and non-Abelian 
quantum Hall states. 



F. The stability of topological orders 

We have analyzed simple topological orders of spinful 
particles using an SU(2) gauge theory. In several occa- 
sions we first considered Laughlin states with conserved 
spin, and then naively generalized the analysis to the 
Rashba spin-orbit coupling that does not conserve spin. 
However, real materials do not feature the SU(2) gauge 
structure and contain perturbations that violate all sym- 
metries of the ideal spin-orbit couplings. We would like to 
explain here using duality why none of these violations of 
spin-conservation or the SU(2) symmetry jeopardize the 
conclusions so far, or why the spin-related topological 
orders are stable. 

We found in the section lHI Al that generic U(l)xSU(2) 
Laughlin incompressible quantum liquids feature quasi- 
particle excitations with fractional quantum numbers 
given by (I36p . regardless of whether the spin-orbit cou- 
pling conserves spin or not. What matters in the deriva- 
tion of that result are only the conditions (A) and (B) 
stated at the beginning of the section IIII Al Specifically, 
the vortex densities p s in all spin channels must be effec- 
tively frozen. Note that the conditions (A) and (B) make 
no reference to the possible conservation of the S z spin 
projection. If S z fluctuates, then the winding numbers n s 
from (|32|) and (|36l) fluctuate accordingly in a state where 
all p s are fixed. The crucial point is that the dual-spin 
of vortices ip in (TTJ) is conserved, even when the spin of 
particles £ in (T74l is not conserved. This is what allows 
the vortex densities p s in various dual-spin channels to 
have quantized incompressible values. 

The conservation of dual-spin can be made explicit 
by emphasizing its nature in the context of the parti- 
cle Lagrangian ([74"]). A vortex with a definite dual-spin 
projection is a topological defect of a particular particle 
spinor component £ s = |£ s | exp(iip s ) in which the phase 
ip s winds by 2irx integer about the singularity. The con- 
servation of dual-spin corresponds to the conservation of 
vorticity in the particle theory (T74"]) . If the particle La- 
grangian is allowed to have only local terms, then only 
vortex-antivortex pairs and other topologically neutral 
deformations of the fields in any given spin channel can 
cost finite action (energy). Furthermore, no local opera- 
tor written in terms of £ s and £] can create or annihilate 
a single vortex, because it would have to qualitatively 
alter the field configuration arbitrarily far away from the 
vortex core. 

A local theory can allow the non-conservation of vortic- 
ity only if it is a gauge theory with deconfined monopoles 
(typically a compact gauge theory). A gauge field cou- 
pled to the matter field is necessary in order to compen- 
sate the macroscopic energy cost of a vortex configura- 
tion in the matter field alone. Then, creating or annihi- 
lating a vortex at a particular point in space and time 
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requires compensation by a monopole gauge-field config- 
uration. In our case there are 2S + 1 flavors of vortices, 
and we could introduce that many compact gauge fields 
X sjl in the model (l?4l to spoil the conservation of vor- 
ticity. Alternatively, we could make the gauge fields X sfl 
non-compact and combine them in the following kind of 
the Maxwell term: 



s\ 



(76) 



This would prohibit individual monopoles, but allow 
monopole-antimonopole pairs in arbitrary two spin chan- 
nels that correspond to dual-spin flips (at the expense of 
the SU(2) symmetry). Both scenarios require dynami- 
cal gauge fields, which are simply absent from the theory 
(|74|) of particles. At best, the particles may be charged 
and coupled to the conventional U(l) gauge field of elec- 
trodynamics (which has non-compact dynamics) , but the 
spin-orbit SU(2) gauge field is fundamentally static so 
that spin-flavored vorticity must be conserved. Note that 
the choice of the spin projection axis is arbitrary in the 
above argument, so the full SU(2) dual-spin symmetry 
can be partially lowered in ([1} only through the external 
fluxes <&o and Qq. This is at the heart of the robustness 
of topological orders against perturbations such as disor- 
der and spin non-conserving collisions. Perturbations do 
not introduce dynamical gauge fields. Topological order 
can be destroyed only via the Landau-Ginzburg part of 
the Lagrangian by perturbations that can overcome the 
TPs gap and make a conventional state more energeti- 
cally favorable. 

Contrast this to the situation of the dual Lagrangian 
of vortices (JTJ) . It does feature a dynamical gauge field 
Bp, which describes particle fluctuations in all spin states 
and provides the means to allow the non-conservation 
of dual vorticity, that is the physical particle spin. One 
weighted trace of B^ is related to charge and must be non- 
compact due to charge conservation, but certain combi- 
nations of B^ eigenmodes can be compact and implement 
spin-changing events via appropriate modifications of (j5J| . 
An example of such a Maxwell term is: 
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which is written using the continuum notation but ul- 
timately defined on a lattice, using the representation 
B^ = diag(m_s^_s )jLl , . . . ,rnsf3s,^) that was introduced 
in the previous section. Generally, the compact cosine 
terms have to reflect the specific ways in which the spin 
is not conserved, and therefore depend on the spin-orbit 
coupling and various spin non-conserving perturbations. 
We have no means to derive by duality the precise form 
of the Maxwell couplings in these general circumstances. 
Nevertheless, we can understand qualitatively the effect 



of fluctuations that are made possible by the compact 
gauge fields in the vortex Lagrangian (p}. 

A compact gauge field can admit 27r-quantized flux 
lines of arbitrary length without any energy cost, im- 
plying that monopoles (ends of semi-infinite flux tubes) 
are allowed. The ensuing fluctuations of quantized flux 
tubes and monopoles lead to the quantization of lo- 
cal density fluctuations in all charge and spin channels 
that couple to the compact gauge fields (in the spirit 
of the Dirac's charge quantization by the existence of 
monopoles). The meaning of density quantization is that 
an isolated volume of space can contain only integer mul- 
tiples of charge/spin quanta. The quanta can be mo- 
bile, but the suppression of smooth density fluctuations 
produces an incompressible ground state. The specific 
compact components of the dual Maxwell term Cm dic- 
tate which particular vortex density fluctuations become 
quantized in the theory (JT|). For example, (|77p would 
quantize all vortex density fluctuations that carry a def- 
inite dual-spin projection on the z-axis. A different kind 
of spin non-conserving dynamics, translated by duality 
to a different Cm, could quantize other dual-spin projec- 
tions, possibly dependent on position or momentum. 

The Landau-Ginzburg part of the Lagrangian © has 
to determine the incompressible vortex densities that 
shape the topological order while respecting the imposed 
density quantization rules. We cannot predict the out- 
come of this without knowing the precise form of Cm , but 
we can rest assured that only the dual-spin densities of 
vortices can be affected. Consequently, the quantization 
(1341) of p s holds against all forms of spin non-conservation 
at least in the TR-invariant Laughlin states, because the 
resulting average dual-spin density is zero and properly 
quantized. However, spin non-conservation affects the 
quantum numbers of excitations. The fractional quan- 
tum numbers of quasiparticles discussed in the section 
IIII Al can be jeopardized, but fractional vortex excitations 
(typically discussed in the FQHS literature and explained 
in the section ITV Al) are protected by the dual-spin con- 
servation. 



IV. GENERALIZATIONS OF TOPOLOGICAL 
ORDERS 

This section demonstrates how the theory ([1} and its 
dual (|74| can describe a broad range of topological in- 
sulating states. We will first develop the formalism for 
describing hierarchical Abelian FQHS and immediately 
extend it to the general Abelian SU(2) states of fractional 
TIs. This formalism will also provide a natural descrip- 
tion of uncorrelated "integer" quantum (spin) Hall states. 
Then we will discuss non- Abelian states and generaliza- 
tions to arbitrary symmetry groups and representations. 
We will demonstrate how the topological term of (|74l) can 
produce non- Abelian statistics of excitations depending 
on the character of low-energy fluctuations, and devote 
a special attention to new topological orders that could 
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be obtained from the Rashba spin-orbit coupling. 

Many important properties of the ground state are de- 
cided by the dynamical (Landau-Ginzburg) part of the 
Lagrangian. They include particle densities, excitation 
spectra, vortex winding numbers and the corresponding 
fractions of elementary charge or spin that are expressed 
in quantum fluctuations. A class of ground states in ex- 
ternal gauge fields features frustration which is resolved 
by nucleating vortices and "binding" them to particles 
in an incompressible quantum liquid. Instead of trying 
to microscopically understand this "binding" process, we 
qualitatively capture its most important outcomes via 
the topological term. The unique feature of (QJ is that 
its topological term is very general and does not by it- 
self specialize to any concrete topological state of mat- 
ter. Instead, the topological term describes the quantum 
entanglement due to particle-vortex "binding" once the 
conventional dynamical properties of the ground state 
are known (or chosen for classification purposes). The 
deducible manifestations of the many-body quantum en- 
tanglement are the fractional statistics of quasiparticles 
and the ground-state degeneracy on Riemann surfaces 
such as torus. Therefore, the topological term alone is 
responsible for describing topological orders, and can be 
used to classify them. We will here discuss only topologi- 
cal orders without asking what dynamical conditions are 
necessary for stabilizing them. 

The section MI Al characterizes some entanglement ef- 
fects in a particular set of incompressible topological 
states whose vortex densities are given by ([Ml) . We can 
now imagine more complicated incompressible states that 
feature multiple low energy excitations labeled by some 
"flavor" quantum number i — 1, . . . , n that corresponds 
to emergent symmetries. The appearance of emergent 
symmetries in the low-energy dynamics is the only kind 
of ground state reorganization that overcomes the very 
restrictive condition without leading to a conven- 
tional state such as superconductor or Mott insulator. 
The resulting possible topological states of matter form 
a hierarchy based on the emergent symmetries, and are 
generally obtained form ([T]) when the background flux $o 
and density &q matrices do not commute. 

More generally, the low-energy dynamics can lead to 
non-Abelian topological orders. Only the topological 
term in ([1} can capture them properly, while the Landau- 
Ginzburg part written there is specialized for the dynam- 
ics that favors Abelian statistics. We will not attempt in 
this paper to generalize this Landau-Ginzburg part be- 
cause its form is known only to the extent allowed by 
the duality transformation of the physical particle theory 
(1741) , and the duality mapping is currently available only 
for the cases with Abelian statistics. Instead, we will ex- 
plore the non-Abelian topological states starting from the 
particle Lagrangian (174")) whose Landau-Ginzburg part 
can be readily constructed by understanding the micro- 
scopic system of interest. The goal we pursue in this 
paper is very modest and limited to the construction 
of effective theories of a few highly-entangled topologi- 



cal liquids, should they be stable ground states. Many 
important issues will be left untouched, most notably 
the characterization of any observable properties. A sys- 
tematic study of fractional non-Abelian states and their 
properties is left for future work. 

The main message will be that the topological spinor 
Lagrangian provides a formalism of much greater flexi- 
bility than the standard CS theory. Its topological term 
can be reduced to an effective CS form by coarse-graining 
in an incompressible quantum liquid, but the CS gauge 
fields may be arbitrary functions of any number of in- 
dependent parameters that are restricted by the sys- 
tem's fundamental symmetry group and its representa- 
tion. These parameters need not be sufficient to generate 
fully unconstrained fluctuations of the CS gauge fields. 
They instead represent the physical low-energy degrees of 
freedom, while the CS gauge fields are merely the math- 
ematical tool that endows the physical excitations with 
non-trivial statistics via the topological term. We will 
show that such a flexibility is necessary to describe new 
topological orders that could arise in the Rashba spin- 
orbit-coupled TIs. 



A. Hierarchical states 

First, we will consider Abelian fractional quantum Hall 
states of electrons in external magnetic fields (without a 
spin-orbit coupling), and obtain their most general CS 
descriptions from (|TJ) . The crucial assumption we need 
to make is that the low energy dynamics exhibits an 
emergent U(l) n symmetry. The number of low energy 
modes, and the level of the effective CS theory, is equal 
to the dimension of the U(I) n symmetry representation. 
We will consider only the minimal n-dimensional rep- 
resentations (conserved quantum numbers), so the vec- 
tor ip — (ip , ip n ) in ([1} will group together n vor- 
tex flavors. The corresponding flavors of particle modes 
are captured by the n-dimensional gauge field matrix 
that depends on n fluctuating eigenmodes. These can 
be either microscopic particles in different internal states 
or any emergent low-energy modes that carry conven- 
tional quantum numbers (charge, spin, band/orbital in- 
dex, etc.). The currents f of particle modes can be 

related to vortex currents j* via a linear transforma- 
tion: 

4 = E ylJ ^ A ^vA (78) 

i 

3 

where the coefficients Y lJ have a dynamical origin in £lg 
and are related to vortex charges along the lines hinted 
in the section Ull El We will also consider a generic static 
U(l) n gauge field coupled to vortex flavors in the 
topological term. This extended gauge field generally 
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depends on the physical U(l) gauge field a M . All space- 
time components of A^ commute with one another, but 
need not commute with B^, in which case we have a non- 
trivial linear relationship ([75)1 . 

By the assumed symmetry, we can choose to work 
in the representation that simultaneously diagonal- 
izes all components of An and $, t . Let us define 



I p 1 exp(i0 4 ) in this representation, and A^ 



di&g(A* , . . . , A"), $ = diagO 1 , . . . , 0"). If the fluctua- 
tions of all densities p l are suppressed, then we can re- 
peat the analysis from the section UlI Al and relate the lo- 
cal measurements of microscopic particle quantum num- 
bers (charges) to the quantized windings of 8i in the low- 
energy vortex field configurations. From (|78|l we find: 



Q' 



3 
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It is convenient to define Z %3 = lisY^ty 1 p> and switch 
to matrix notation Q — Zm, where the vectors Q = 
(Q , . . . , Q n ) and m = (m , . . . , m n ) have integer com- 
ponents. The detection of a physical particle excitation 
in the flavor state k corresponds to the vector Q{k) with 
components Q l = di.k, and is related to a particular 
combination of dual-vortex winding numbers contained 
in the vector m(k). We find m(k) = Z^ 1 Q(k), so that 
all matrix elements of Z~ l must be integers. A particular 
compliant matrix Z characterizes the topological insulat- 
ing ground state and determines the fractional quantum 
numbers SQ = (SQ 1 , . . . , 8Q n ) of various quasiparticle 
excitations enumerated by integer- valued vectors I: 



SQ = Zl 



{z- 1 )^ , r e z 



(80) 



The constraints on Z still leave significant freedom 
for the matrix Y. The main restriction on Y comes 
from the definition Z ij = 2ttY v ' (fp p 3 which implies 
that Z ij /Y 13 can depend only on the index j, that is 
Y l i = Y lj (Z 1 ^ /Z 1 !). Otherwise, we are free to impose 
other requirements without any loss of generality. For ex- 
ample, we can require that the matrix Y be orthogonal. 
This amounts to a choice of normalization for tjj 1 . The 
ground state is characterized by incompressible vortex 
densities: 



1 Z 1 ? 
2tt^ Y^ 



(81) 



By fixing p % , the currents (1751 and the topological term 
in the Lagrangian (J]) reduce to: 



(82) 



This can be expressed in the matrix form by introducing 
a diagonal matrix R = diag(^ 1 p 1 , . . . , 4> n p n ) and vectors 



6^ = (^ 1 ,...,W,JPM = fc.-.,J?J: 



±-Ze» vX dM 
2ir 



(83) 



C t = -e^ x b T .Rd„b, -te^ x A^Rd„b, 



The currents of particle modes can be directly repre- 
sented as fluxes of a different set of CS gauge fields 
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(84) 



2ir 



(79) where = Zb lll and the matrix K and vector q arc 



defined by: 
K = 2tt{Z' 



- X \ T RZ~ X 



a til 



T = 2irA T RZ~ 1 . (85) 



The topological Lagrangian in (1841 now has the standard 
CS form, with the coupling matrix K being symmetric 
by definition. Further restrictions on K and q follow 
from the requirement that j PM be the physical particle 
modes with conventional quantum numbers, which can 
be created or annihilated by local combinations of physi- 
cal electron operators. It can be easily seen from (|75|) and 
subsequent definitions that measurable (integer) quan- 
tum numbers of physical particle modes correspond to 
flux quanta of the c^, gauge fields. One of the particle 
modes directly corresponds to electrons and thus carries 
the U(l) charge that couples to the external gauge field 
a M and has fermionic statistics. Let this mode be la- 
beled by i = 1 in our representation, so that the vec- 
tor q is given by q — (1, 0, 0, . . . , 0). The other modes 
are normally neutral (particle-hole) and hence must have 
bosonic statistics. The CS self-coupling implements the 
exchange statistics of excitations via the matrix K, and 
since the physical modes are flux quanta of , the matrix 
elements of K must be integers. Specifically, K 11 must 
be odd and K 11 , i > 1 must be even. Other choices for q 
and K are appropriate for systems with different types of 
particle modes. Note that (Z -1 )'^ £ Z is not enough to 
make K tJ £ Z, the latter requirement further constrains 
the possible vortex densities p 1 in quantum Hall states. 
From the equation of motion one obtains the filling factor 
v = q T K~ 1 q. 

The literature on fractional quantum Hall effects^ 
generally considers a different kind of fractional excita- 
tions than the ones captured by (|8"0"|) . These excitations 
are "fractional vortices" whose currents — (J 1 , . . . , j") 
minimally couple to the CS gauge fields c^. Their U(l) 
charge Sq and statistical angle 89 are given by: 



Sq = q T K- 1 l 



(86) 



where I is any vector of integers in the flavor space. The 
quasiparticles of (I8TH) minimally couple to the gauge field 



A^ and thus are dual to vortices. 
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It should be emphasized that the topological La- 
grangian C t and the CS theory as its special limit do 
not by themselves determine the character of fluctuations 
that are detected in any particular experiment. The na- 
ture of low-energy modes is (also) shaped by the Landau- 
Ginzburg part of the Lagrangian. The eigenstates of 
the appropriate many-body Hamiltonian always carry an 
integer-quantized total charge, which however may be 
spatially distributed in fractional lumps. The topological 
Lagrangian only imposes some constraints on what kinds 
of lumps are possible and how their relative locations af- 
fect the many-body wavefunction's phase. The dynamics 
of lumps is beyond the topological term's reach, but af- 
fects the statistics of measurement outcomes. In that 
sense, one cannot easily make predictions about what 
kind of excitations would a particular experiment be sen- 
sitive to, fractionalized particles, vortices, or some other. 
A prediction we can make is that if an experimentalist 
successfully localizes a single fractional vortex, the ob- 
servable amount of charge in its vicinity will be given by 



The effective Lagrangian in (f84|) is the most general 
CS theory of hierarchical Abelian quantum Hall states. 
We have seen that such topological states can resolve 
the frustration of electron's kinetic energy at virtually 
any fractional particle density of v particles per flux 
quantum, provided that the low-energy dynamics spon- 
taneously develops multiple internal degrees of freedom 
for quasiparticle and vortex excitations. Analogous hi- 
erarchy of fractional states can now be constructed for 
TR-invariant topological insulators. The Abelian hierar- 
chical states of spin 5 particles feature the SU(2)™ sym- 
metry, with the (25 + l)-dimensional representation of 
the SU(2) subgroup. A (25 + l)-component sub-spinor 
= (ip l '~ s , . . . ,if>' l ' S ) is needed for each vortex flavor i 
to determine vortex "charge" 2\a an d "spin" currents 
analogous to (|26p: 



Jvfi 



J 1 ' 



(87) 



We can allow a flavor-dependent SU(2) flux $q to ensure 
the proper TR transformations. Then, we can generalize 
the physical particle charge and spin currents (|T8[) as: 



J 



vX B ( Y lj i j 

u v \ l 00 J V X 

vX f) () " I ' 
°u I I a0 JvX 



pa, 

■yij Jjb 
1 ab "vA 



We used here the Einstein's notation for all indices. The 
spin-orbit coefficients Y^ a and Y^ must be either zero, or 
operators that change sign under TR. If in addition to the 
emergent SU(2) n symmetry the low-energy dynamics fea- 
tures mutually commuting emergent gauge fields A 1 ^ , and 
their fluxes <I>^ in all flavors, then the topological term C t 
has an even higher U(1) < - 2S+1 '" symmetry and reduces to 



one of the CS theories discussed in Ref42 when the den- 
sity fluctuations are small. It is also straight-forward to 
construct the CS theories of hierarchical Abelian frac- 
tional quantum Hall states for electrons that experience 
both a strong magnetic field and spin-orbit coupling. 

The TR-invariant topological gauge theory of spin 5 = 
i particles obtained in this manner can be written in the 
BF form: 



J z — 
A = 



1 



(89) 



An 



4tt 



T jz 



The superscripts c and s label the CS gauge field vectors 
that represent charge and spin currents respectively. 



B. Non- Abelian states 

Low-energy dynamics can support conditions for non- 
Abelian fractional incompressible quantum liquids. It is 
not presently clear how to describe such conditions in 
the vortex Lagrangian (fT]), so we will construct the non- 
Abelian effective theories in the language of the particle 
Lagrangian (|74|) . The prominent fluctuations of the n- 
component particle spinor field £ could be captured in 
certain topological states by m fluctuating phases ip a , 
where a £ {1, . . . , m} and m < n. A fairly general form 
of such fluctuations can be written as: 



£ = e l(p *> 



( f 1 



(90) 



in some representation, where r] a are a set of linearly 
independent Hermitian n x n matrices, and f a are non- 
fluctuating complex or Grassmann amplitudes. We will 
be interested in non-commuting r) a (the above low-energy 
fluctuations generated by mutually commuting rf could 
produce only Abelian topological orders). If we define 
the gauge field matrix 

= -z(V W V^ a (91) 



o 

we can write: 

= iZrf . (92) 
The charge and spin currents (|75p are: 

j„ = tHz, - AJZ = tr[(Z„ - A^)] (93) 
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In the limit of suppressed amplitude f a fluctuations, the 
topological term from (j?4")) becomes: 



•Cpt — 



-tr 



up to a total derivative, where 



e^\d v Z x -iZ v Z x ) 



(94) 



(95) 



It was pointed out in the section Ull El that £ pt is gauge- 
invariant, but written in the natural gauge where the 
external gauge field _4 M is directly derived from the fixed 
form of the Rashba (or other) spin-orbit coupling. Mak- 
ing (|94p manifestly gauge-invariant requires replacing Z^ 
by Z M — V M , where V M is a static flux- less background 
gauge field coupled to vortices. The dependence of 



/ I/ 1 ! 



t _ Jv a v" 



-lip r\ 



(96) 



n\2 
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on ip a generally seeps into the Lagrangian (|94l) . which 
therefore is not a pure gauge theory. Nevertheless, the 
kinematics shaped by features excitations with non- 
Abelian fractional statistics. 

The fluctuations of tp a contain both singular and non- 
singular components, and the topological kinematics of 
the singular ones is produced entirely via the non-Abelian 
gauge field Z„. Being interested only in the qualitative 
aspects of topological kinematics, we are tempted to con- 
vert (|94[) into an approximate pure gauge theory. This 
can be done by coarse-graining. We will integrate out 
the short length-scale fluctuations in the path-integral 
that average out the non-singular fluctuations of ip a and 
allow the CS flux to diffuse. We need to identify the in- 
dependent degrees of freedom <p a that adequately span 
the prominent low-energy fluctuations, and this amounts 
to choosing the proper set of matrices rj a and amplitudes 
f a that match the dynamics (we are not prescribing the 
method to do this here). Then, we can formally treat 
ip a as independent random variables. We will also re- 
place ££t in (|M|) by its average under the assump- 
tion that fluctuations are abundant and chaotic at short 
length and time scales. This amounts to neglecting the 
subtle microscopic correlations between ££t and Z^ (that 
is (p a and d^ip a ). 

We will consider two examples of the above coarse- 
graining procedure. The first example is a level-2 hierar- 
chical state of electrons in magnetic field: 



f 1 
f 



(97) 



where <p — (cp x ,cp y ,0) and er is the vector of Pauli matri- 
ces <j a . This two-component spinor represents prominent 
fluctuations of two particle modes in some unspecified 



quantum liquid. The gauge field (|92l is: 



1 



sin(2|</?| 



2\<p 
sin(2|y> 



m 



(98) 



sin \ip\ 



(if x d^ip)o- 



= Z^,d^)a a . 

All components Z" are non-zero for generic values of 
ip, so the gauge field is "microscopically" non-Abelian. 
Integrating out the short-wavelength fluctuations of ip 
will coarse-grain the gauge field in a manner that is 
hard to extract analytically due to the non-linearity of 
Z®(ip, d^ip). However, we can rest assured that even the 
coarse-grained gauge field is non-Abelian because it is 
traceless while being dependent on two angles per point 
in space-time. One of these two angles would have to 
control the fluctuations of tr(Zp) in order for Z^ to be 
Abelian in its two-dimensional representation. 

We can approximately coarse-grain £ and ££t by aver- 
aging them over the independent local random variables 
f x (r) and fy{ Y )- Assuming that the probability distri- 
bution is symmetric under ipi — > — ipi, we find: 

(cos \<p\) ( f f2 



(0 



(99) 



and 



tv _ \ ff±\ff , I/ 1 ! 2 - l/ 2 l 
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m = 



(cos(2|( y a|))cr z (100) 



2 2 
+Re (ff 1 *) (cos 2 MK + Im (ff 1 *) (cos 2 \<p\)^ . 

This can be now substituted in (|94|) to obtain a pure 
non-Abelian effective CS theory. The amplitudes /' can 
be determined in the spirit of the section Ull A[ pending 
some information about the microscopic physical excita- 
tions that ([97]) represents. For illustration purposes, let 
us make the simplest assumption that the isolated phase 
windings in ip x or ip y correspond to microscopic (non- 
fractional) bosonic vortices. If we generate only one type 
of vortices, via ip x for example, then the other phase 
ip y does not have singularities and does not contribute 
to the topological term. The gauge field created 
solely by (p x becomes Abelian, Z^ — a x d^Lp x = C,^a x . In 
the present representation, this Abelian Z^ and the most 
general Oq" 1 = t + \Tga a produce the following effective 
theory from ([94]) : 

W vX W x ) (101) 
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We can now deduce the constraints on f % by reading out 
the statistics of excitations from C x t and the analogous 



to(|/ 1 | a + |/ 2 | 2 )+IS , Re(/ a / 1 *) (cos 2 |^|) 
t (\f\ 2 + \f\ 2 )+T»Im(ff*){co S 2 \<p\) 
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The (even) integers m x and m y characterize the topolog- 
ically ordered ground state whose fractional vortex exci- 
tations have mutually non-Abelian statistics. 

It should be noted that the non-Abelian CS theory 
(|94|) applied to the present example (|97|) is not a true 
gauge theory even though it looks like one. The actual 
fluctuations are generated by two local variables rather 
than a real SU(2) gauge field that requires three local 
variables (for three "gauge boson" modes) to fully explore 
its phase-space. The measurable currents associated with 
the two particle modes are obtained by substituting (|97[) 
and (J55JI into ([55]). 

A variation of this example is the non-Abelian incom- 
pressible quantum liquid of spin S = h particles whose 
spinor components in (|97p are identified with the spin S z 
amplitudes, f 1 — fa and f 2 — /j.. This state is generally 
not TR-invariant, but it can be TR-invariant if certain 
conditions are met. The necessary conditions for the TR 
symmetry are that the averages (cos(n\(p\)) vanish, and 
one of the amplitudes fa or f± be zero. The first condi- 
tion is naturally obtained from the wave-like oscillations 
(pi ~ pr — tot. Then we would get (£) = and oc 1 
consistent with the TR symmetry. 

The second example is a TR-invariant topological in- 
sulator of S = 1 particles such as p-wave Cooper pairs. 
We can construct a symmetric non-Abelian pure gauge 
theory starting from the normalized spinor: 




(102) 



in the representation that diagonalizes the external spin- 
orbit SU(2) flux $o- The matrices j a are here the three- 
dimensional representations of the three SU(2) genera- 
tors. The fluctuations generated by ip a span the full 
three-dimensional spin space in a manner that treats dif- 
ferently the xy plane and the z-direction, but respects the 
TR and xy rotation symmetries. Therefore, the symme- 
try of fluctuations is consistent with the spin U(l) sym- 
metry of the model, which is anyway established by the 
presence of the external spin-orbit SU(2) flux <£> = &qJ z ■ 
By following the same procedure as before, the coarse- 
grained gauge field and density matrices are found to be: 

Z M = Z>*,^,^) 7 a , (&)=a+Ptf , (103) 

where a and f3 are numerical constants (note that (j z ) 2 ^ 
1 in the S = 1 representation) . No symmetry is violated, 
and one obtains a pure non-Abelian gauge theory by sub- 
stituting this into (|94l) : 



la 

£ pt = rtr 



(104) 



tr 



{2„,e^}{*£,*8} 



We will not pursue in this and subsequent examples the 
quantization of constants such as a and /3 that repro- 
duces the correct statistics of microscopic excitations. An 



n-dimensional gauge field representation can be Abelian 
only if it depends on up to n parameters and one them 
controls the tr(i^). But, the gauge field here belongs to 
a three-dimensional representation and depends on three 
parameters while being strictly traceless, so it is guar- 
antied to remain non-Abelian after coarse-graining. Its 
scalar components can be regarded as independently 
fluctuating fields after coarse-graining, so (|104[) is a true 
SU(2) non-Abelian CS gauge theory, which respects the 
TR symmetry 

We would like to remark here that the obtained non- 
Abelian gauge theories have a self-coupling of the dynam- 
ical gauge field that differs from the usually discussed 
(TR-violating) CS form of the Yang-Mills gauge field: 
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(105) 



This is gauge-invariant up to the boundary term, and the 
Wess-Zumino-Witten ter m 65 i 66 whose gauge-invariance 
requires the quantization of the constant k. In contrast, 
the form (I94|) is gauge-invariant because it contains a 
(hidden) background gauge field V M bundled into Z^ — V M . 
Even though we could set the gauge V M = 0, the funda- 
mental presence of V M is very much in the spirit of the 
duality between ([1]) and (f74|) . The formally analogous 
topological terms of ((TJ and (jT4"]l are essential for repro- 
ducing the correct drift currents via equations of motion. 
The CS coupling (jl05[) might require a different kind of a 
topological term in the spinor Lagrangian, possibly tai- 
lored to a different kind of physics than discussed in this 
paper. 



C. Non-Abelian topological orders from the 
Rashba spin-orbit coupling 

The Rashba spin-orbit coupling H so = vz(p x S) in 
quantum wells does not conserve any fixed spin projec- 
tion, but respects a special dynamical symmetry that we 
discussed in the section IIII CI If no perturbations or in- 
teractions spoil this symmetry, it will define the quantum 
numbers of excitations. Our goal here is to rudimentar- 
ily explore the highly entangled incompressible quantum 
liquids consistent with this symmetry, and establish that 
they host quasiparticles with a non-Abelian fractional 
statistics. The latter indicates a non-trivial topological 
order, which is automatically robust against sufficiently 
weak perturbations even if they spoil the symmetries of 
the ideal Rashba-coupled Hamiltonian. 

Let us focus on spin S = \ particles. If S z were con- 
served, we would expect the spinor configurations 



f(r) - ( Cte ^ (r) 



(106) 



to be relevant at low energies when the density fluctu- 
ations are small in comparison to the average densities 
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|Cf | 2 and |CjJ 2 . This applies to quantum Hall and spin- 
Hall liquids, and leads to our earlier constructions of the 
Abelian CS theories, albeit in the language of vortices 
rather than particles. The S z conservation here defines 
the two eigenmodes that carry the opposite S z spin pro- 
jections, and the independent fluctuations of these modes 
at fixed densities are captured by the above spinor con- 
figurations. 

In contrast, the translation symmetry and the dynam- 
ical symmetry of the Rashba spin-orbit coupling create 
particle modes that carry a conserved momentum p and 
a conserved "helical" spin projection (z x p)S. Generic 
excitations are superpositions of these modes that we can 
describe by the following spinor field configurations: 



e±(p) 



/±(p) 

V2 V ±e 



/±(p)e 



( 10? ) 

where (f> p = arg(— p y + ip x ) is the angle of the in-plane 
direction perpendicular to p according to the right-hand 
rule, and er is the vector of Pauli matrices. The subscript 
± indicates the sign of the helical spin projection. The 
real-space field configurations £(r) 

e ± w = J (02 eipr ^±(p) ( 108 ) 

are characterized by the densities 



/±(P)/±(P') ■ 



(2tt) 2 (2tt) 2 
1 + g^p-'V) 



An incompressible quantum liquid that conserves the he- 
lical spin projection will have constant uniform densities 
|£±( r )| 2 = |C±| 2 in both ± sectors, so that 



d 2 q 1 



a»Wq-0q-p) 



(27T) 



-/±(q)/i(q-p) = |27rC±| 2 5(p) . 

( 11Q ) 

This represents one constraint on each of the two complex 
amplitudes f± per point in space-time. Therefore, the 
allowed configurations of £±(r) arc determined by one 
real function ip±(j) in each ± sector. The complexity 
of (|110l) will unfortunately prevent us from identifying 
the independent variables <p±(j), but we will still gain 
some insight about the quantum states shaped by their 
fluctuations. 

The TR-invariance and other symmetries impose no 
restrictions on the incompressible densities C + and C_. 
The field configurations £+ and £_ are generally not de- 
generate, and the dynamics could freeze the fluctuations 
of one of them. Let us consider such a special case, where 
for example only £+ is important. Having a fixed density 
in an incompressible quantum liquid allows us to 
generate the low-energy fluctuations by applying a uni- 
tary operator to a fixed spinor: 



£+(r) = l> + [r;^+(r)] 







(111) 



The CS gauge field (J92]) is given by Z M = -i{dJj + )U\ 
as a function of the independently fluctuating variables 
ip+(r). All aspects of topological order in an incompress- 
ible quantum liquid arc then captured by the effective CS 
Lagrangian (|94|) expressed in terms of this gauge field, 
and by the knowledge of low-energy particle modes that 
restrict the gauge field configurations. The effective the- 
ory also contains , but the TR-invariance ensures that 
it can be replaced by = 5IC+I 2 x 1 after coarse- 

graining. The winding numbers of vortex excitations de- 
termine the quantization of |C+| 2 that appears in the CS 
theory, and hence the fractional amounts of charge and 
spin. The ensuing topological orders can be classified by 
a Laughlin-like sequence of states and hierarchical con- 
structions. However, the statistics of quasiparticles and 
the ground state degeneracy on Riemann surfaces are ex- 
pected to differ from the ordinary Laughlin and hierar- 
chical quantum Hall states. The statistics of excitations 
can be found by considering the adiabatic evolution of a 
quasiparticle's spinor as it moves on a closed path around 
another quasiparticle. It depends on the CS gauge field 
in the vicinity of a vortex, which in the present case 
gives rise to a non- Abelian Aharonov-Bohm effect and 
a non- Abelian fractional statistics once |C+| 2 is fraction- 
ally quantized. Unfortunately, we cannot reach more de- 
tailed conclusions at this time. The characterization of 
non- Abelian topological orders is an open problem, since 
the steps outlined here cannot be carried out without a 
better understanding of the non- Abelian duality. 

Nevertheless, we can at least get a glimpse of a £ + 
vortex and appreciate why the quantum liquids of such 
vortices produce TR-invariant non-Abelian topological 
orders. We do not know how £7+(r) in (|111|) depends 
on (p + (r), but we can establish that it is locally a lin- 
ear combination of the unit-matrix, a x and a y , which 
does not include a z . This follows from the fact that the 
Pauli matrices are introduced in (|107|) and (|108l) only via 
the products per with a two-dimensional vector p. Any 
unitary operator of this kind can be written as 



U 4 



(112) 



In our case, fo( r ), <Ar( r ) an d <p y (r) are mutually- 
dependent functions derived non-locally from a single in- 
dependent function ip + (r). After isolating out the trivial 
charge fluctuations tpa, the rest of this operator is the 
same as the one in (|9~7| , so the traceless part of the cor- 
responding CS gauge field = iZ^ + takes the form 
in terms of the vector cp = (ip x , f y ,0): 



sin(2M) 



2M 



sin(2M) 
2M 



sin \tp\ 



(113) 



{ifi X d^)(T 



An elementary £ + excitation that carries momentum 
p has an amplitude /+(p') oc 5(p' — p), so its field con- 
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figuration is: 



£i(r) = 



V2\ 



,i<t> P 



(114) 



This mode carries both a charge and spin current. 
Charge currents must be absent in a TR-invariant ground 
state, so let us consider a superposition of two elementary 
£ + modes that carry counter-propagating currents: 



= - 



1 



1 



-ipr 



oi(pr)[(zxp)<r] 



(115) 



The corresponding CS gauge field d^2 — is: 

Z, L = [(z x p)cr]pi , (116) 
and it can be reproduced by 



<A)=0 , y = (pr)(zxp) 



(117) 



from the expression (|113[) . It can be seen either from (|75|) 
or ()93j) that the charge and spin current densities carried 
by this excitation are: 



3i =0 



(118) 



meaning that only the helical spin projection perpendic- 
ular to the momentum p is carried in the direction of p. 
The charge and spin densities are: 



Jo 



= &a = 1 



(119) 



JS* a = = i[zcos(2pr) 



p sin(2pr) 



The number (charge) density is uniform and thus con- 
sistent with incompressibility. The spin density appears 
non-uniform, but it should be noted that it oscillates in 
both space and time. We have been emphasizing only 
the spatial dependence, but pr in fact stands for pr — uit, 
where the frequency ui is by magnitude equal or greater 
than the excitation gap A in the incompressible liquid 
state. These oscillations are removed by coarse-graining 
over time intervals larger than St ~ A -1 . The ground 
state formed by the zero-point quantum fluctuations of 
many modes such as (| 1 1 5[) will not have any spin texture 
or charge currents. 

However, as long as we are allowed to superpose only 
the £+ modes, we cannot eliminate the spin currents. 
These currents are allowed by the TR symmetry, but 
must not have an open-ended flow in equilibrium. There- 
fore, a stable ground state of fluctuating £ + must in fact 
feature spin current loops. Any static arrangement of 
such loops is a TR-invariant vortex lattice, which could 
be stabilized if the particles were bosonic. The ap- 
propriate configurations of Z^ in such a vortex lattice 
would have to involve a vector field (p(r) that accord- 
ing to (|117p changes along any vortex-centered loop as 



ip(r) ~ (pr)(z x p), where p is locally tangential to the 
loop and quantized to make <p(r) single- valued. This in- 
evitably generates a quantized u z flux tube at the vortex 
center via the {up x d^ip)cr part of (jTT3J) . The presence 
of an external SU(2) flux <I>o oc a z can naturally generate 
this type of vortex configurations. A topologically or- 
dered and highly entangled SU(2) quantum liquid state 
is obtained upon the quantum melting of this vortex lat- 
tice. 

This qualitative picture hints that the duality transfor- 
mation involving non-Abelian fields has a lot in common 
with the Abelian case. However, an elementary vortex 
cannot be captured by a purely Abelian Z M ~ a z config- 
uration according to (|113[) , so it cannot have an Abelian 
statistics. As an example, consider a single vortex with 
concentric circular spin supercurrents. The spatial Zi 
component taken in the direction perpendicular to the 
local spin-current flow p is bound to produce di<p paral- 
lel to <p (because tp always points in the radial direction) 
and thus Zi oc tper. On the other hand, the component 
Zj in the direction of p is a linear combination of all a a in 
the presence of the current flow curvatures (which make 
djip not parallel to ip) and hence does not commute with 
Zi. One might say that a dynamical SU(2) gauge field 
can emit its own "gauge boson" quanta and provide a 
flux-feedback to the externally imposed flux. 

The kind of a TR-invariant topological order explored 
here is different than any one that can be described us- 
ing the standard CS theories with proper gauge fields. 
Clearly, the types of many-body entanglement discussed 
in this section can be extended to higher spin represen- 
tations, hierarchical states, inhomogeneous quantum flu- 
ids, etc. Fluctuations can be generated by virtually any 
unitary operator that depends on distributed parame- 
ters and acts on the n-component spinors of fixed am- 
plitudes. Each parameter can be distributed over any 
complete set of generalized coordinates and coupled to 
any coordinate-dependent generator of transformations. 
The Lagrangians |T]) and (|7"4"]) provide the formalism to 
describe all of these possibilities for topological states of 
matter, at least in principle. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we constructed a topological field the- 
ory of particles in two spatial dimensions whose charge 
and spin couple to the external electromagnetic and spin- 
orbit fields. This theory provides a universal description 
of both conventional and topological phases, being an 
extension of the Landau- Ginzburg Lagrangian that im- 
plements a state-dependent quantization of classical dy- 
namics via a topological term. The added topological 
term is innocuous in phases where particles are well lo- 
calized (e.g. Mott insulators) or very mobile (supercon- 
ductors and Fermi liquids). However, quantum insula- 
tors in which the external magnetic field or spin-orbit 
coupling frustrate particle motion are qualitatively af- 
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fected by the topological term and acquire quasiparti- 
cle excitations with fractional quantum numbers and ex- 
change statistics. These incompressible quantum liquids 
are generalized quantum Hall states with topological or- 
der (many-body quantum entanglement extending over 
large distances). 

The topological field theory was written in two mu- 
tually dual and physically equivalent forms: ([73)1 de- 
scribes the physical particles directly, while ([IJ describes 
vortices, the topological defects of particle field config- 
urations. The particle Lagrangian (|74l) must transpar- 
ently capture the known dynamics of particles in external 
magnetic and spin-orbit fields, but cannot alter their ex- 
change statistics in smooth field configurations that dom- 
inate the path integral of conventional phases. Therefore, 
any statistics-altering topological term may couple only 
to topological defects. Unfortunately, no experimental 
evidence of vortex dynamics in the materials of interest 
is available to guide the construction of such topological 
terms. Instead, we had to resort to the dual theory of vor- 
tices (JJJ and construct its topological term first. Topolog- 
ical defects of the dual theory correspond to physical par- 
ticles whose dynamics is measurable. We constructed the 
dual topological term by requiring that the dual action 
be stationary when its field configurations reproduce the 
classical equations of motion for constant drift currents 
of particles in external magnetic and spin-orbit fields. 
The appropriate equations of motion were deduced from 
the generic model Hamiltonian of two-dimensional topo- 
logical insulators that exhibit the quantum spin-Hall ef- 
fect. The Rashba spin-orbit-coupled electrons in realistic 
quantum wells are modeled by exactly the same type of 
Hamiltonian, with the same symmetries and qualitative 
structure. Therefore, by using symmetries we achieved 
an indirect experimental justification of the topological 
field theory. The full picture of particle dynamics and 
topological kinematics is obtained only by establishing 
the field-theoretical duality between the particle (1741) and 
vortex (Q]) Lagrangians, which constrains their relative 
form. 

The vortex Lagrangian ([1} is actually a generalization 
of (dynamically enhanced) Chern-Simons gauge theories. 
Whenever the dynamics of a quantum state can be ap- 
proximately captured by an XY model derived from ([1}, 
its topological properties are similarly captured by an ef- 
fective theory derived from (JXJ) that has a Chern-Simons 
structure. A Landau-Ginzburg theory is already an ab- 
straction that gives up the "irrelevant" microscopic de- 
tails of the system in order to focus on its universal prop- 
erties based on symmetries. The true usefulness of a 
Chern-Simons theory is achieved when even the symme- 
tries are regarded as "irrelevant" in order to focus on 
topological orders. A Chern-Simons theory can be writ- 
ten only after assuming some symmetry, such as U(l)", 
but the topological order that it captures is robust even 
when perturbations remove this symmetry. The topo- 
logical term of the Lagrangian (JJJ) is more versatile. It 
does not assume any emergent symmetry, but rather lets 



the Landau-Ginzburg part determine the symmetries and 
densities in the ground state, which in turn dictate the 
fractional statistics of quasiparticle excitations. 

This paper accomplished several goals. First, we ex- 
plained the construction of the new topological field the- 
ory and its dual. Several known results from the litera- 
ture were re-derived and adapted to the specific features 
of the present theory in order to show that it is consistent 
with the dynamics of systems we wish to model, and that 
it adequately generalizes the Chern-Simons theory in the 
cases of standard U(l) quantum Hall liquids. 

Second, we introduced a formalism based on the SU(2) 
gauge symmetry that can describe any spin-orbit cou- 
pling and view it as the origin of a generalized quantum 
Hall effect. The SU(2) "Hall effect" has a richer phe- 
nomenology that its U(l) counterpart and need not lead 
to the quantization of spin-Hall conductivity. In partic- 
ular, we showed that the Rashba spin-orbit coupling has 
a dynamical symmetry that differs from the conventional 
U(l) symmetry of the standard quantum Hall effect, but 
nevertheless creates fractional ground states with topo- 
logical degeneracy on a torus when interactions stabilize 
an incompressible quantum liquid. The simplest such 
topological orders are of the Laughlin kind and do not 
depend on whether S z is conserved or not. We further an- 
alyzed the quantum numbers of quasiparticle excitations, 
and identified the nature of charge and spin fractional- 
ization in relation to the symmetries of the ground state. 
We did this for any combination of external magnetic and 
spin-orbit fields acting on particles with arbitrary spin. 

Third, we demonstrated the ability of the proposed 
topological field theory to handle a broad spectrum of in- 
compressible quantum liquids with distinct robust topo- 
logical orders. We elucidated the construction of hier- 
archical incompressible liquids, mostly by focusing on 
the Abelian quantum Hall states classified by the U(l) 
Chern-Simons theories. The analogous SU(2) descrip- 
tions of hierarchical Abelian spin-Hall states were ob- 
tained in a straight-forward fashion. Then, we consid- 
ered a few examples of low-energy fluctuations that pro- 
duce non-Abelian fractional statistics in incompressible 
quantum liquids. For instance, if such a quantum liquid 
features a triplet of low-energy modes with SU(2) sym- 
metry, its topological order is described by an effective 
SU(2) Chern-Simons gauge theory that can be derived 
from the topological term of (JlJ. In general, the effec- 
tive theory superficially has a Chern-Simons form, but 
its non-Abelian gauge fields can be constrained to a sub- 
space of all possible configurations through their depen- 
dence on a small number of parameters that generate the 
low-energy fluctuations. We specifically constructed an 
effective field theory of a highly entangled non-Abelian 
incompressible quantum liquid that respects the time- 
reversal symmetry and takes advantage of the Rashba 
spin-orbit coupling to lower its energy. 

The proposed topological field theory is much more 
general than the Chern-Simons theory, and its clean 
structure based on representations and symmetries 
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might provide a broad classification scheme of topo- 
logical orders. The present Lagrangian can be ex- 
tended to any representation of arbitrary emergent 
U(l) ni xSU(2)" 2 xSU(3)™ 3 . . . symmetry groups in which 
particles or vortices are coupled to an external (non- 
Abelian) "electromagnetic" flux of arbitrary direction in 
the space spanned by the symmetry generators. The 
symmetry determines the character of topological orders, 
the representation is related to the available low-energy 
particle/vortex modes and the (non) conservation of their 
quantum numbers, while the flux direction determines 
the mutual statistics of modes. There are many possibil- 
ities that will be explored in future work, together with a 
systematic analysis of observable phenomena such as the 
properties of protected boundary states. 

The topological term constructed in this paper could 
be viewed as the lowest-order member of a sequence of 
topological terms that contain higher powers of flux, (co- 
variant) derivatives and field operators. This sequence 
must be restricted by symmetries and its higher order 
members might lead to additional topological orders if 
they are macroscopically relevant. Further generaliza- 
tions to higher spatial dimensions and other types of mat- 
ter fields would involve constructing topological terms 
that reflect all possible topological defects and their 
symmetry-allowed couplings to external gauge fields. 



Finally, the proposed theory also describes the dy- 
namics of particles or vortices via its Landau-Ginzburg 
part. We demonstrated that the topological term is state- 
independent, but the topological order that it creates 
depends on the dynamically stable density of particles 
and vortices. With this level of description, we can ask 
whether the written theory could be used to chart univer- 
sal phase diagrams that contain both conventional and 
topological states of matter. This question is very much 
worth exploring because no method is currently available 
to solve such problems, other than the numerical exact 
diagonalization of systems with a few particles. 
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